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ABSTRACT 

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 /im-selected 
galaxies at z < 0.5 (z = 0.24) from the Hersc/jeZ-ATLAS survey. We complement our Her- 
schel 100-500 /im data with UV-K-band photometry from the Galaxy And Mass Assembly 
(GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic 
SEDs for a representative sample of 250 /im selected galaxies spanning the most recent 5 
Gyr of cosmic history. We derive estimates of physical parameters, including star formation 
rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at 
z < 0.5 has a far-infrared luminosity in the range 10 10 - 10 12 L Q (SFR: 1-50 M Q yr -1 ) 
thus is broadly representative of normal star forming galaxies over this redshift range. We 
show that 250 /im selected galaxies contain a larger mass of dust at a given infra-red lumi- 
nosity or star-formation rate than previous samples selected at 60 /im from IRAS. We derive 
typical SEDs for //-ATLAS galaxies, and show that the emergent SED shape is most sen- 
sitive to specific star-formation rate. The optical-UV SEDs also become more reddened due 
to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red 
templates. They may therefore be most appropriate for inferring total IR luminosities from 
moderate redshift submillimetre selected samples and for inclusion in models of the lower 
redshift submillimetre galaxy populations. 

Key words: Galaxies: starburst 



1 INTRODUCTION 

In the past couple of decades, our understanding of the Universe 
has flourished as a result of our new-found ability to observe in 
almost all regions of the electromagnetic spectrum. This advance 
is in no small part due to our ability to associate observations at 
different wavelengths with particular astrophysical phenomena and 
link them together, making modern astronomy truly pan-chromatic. 
By observing an astronomical source at multiple wavelengths, we 
may piece together its spectral energy distribution (SED), and by 
comparing the observed SED to models, we may infer the physical 
properties of the source (or sample of sources) that we are studying. 

At ultraviolet, optical and near-infrared wavelengths the SED 
of the average galaxy is dominated by emission from stars; there 
are many tens of different models to which we may compare our 
observations, in the hope of understanding the stellar components 
of astrophysical sources (e.g. Bruzual & Chariot, 2003, Fioc & 
Rocca-Volmerange, 1997, Vazquez & Leitherer, 2005, Anders & 
Fritze-von Alvensleben, 2003, Jimenez et al. 1995, 2004, Maras- 
ton, 2005, Pietrinferni et al. 2004). Such SED model analysis may 
be used to determine basic properties of a galaxy's stellar compo- 
nents, such as its age, metallicity or stellar mass (see e.g. Carter et 
al., 2009, Smith & Jarvis, 2007, Collins et al., 2009, Pacifici et al., 
2012, PfoiT, Maraston & Tonini, 2012). 

While the ultraviolet to near-infrared emission tells us about 
the stellar content of a normal galaxy (subject to correcting for at- 
tenuation by dust of the different stellar components; Chariot & 
Fall, 2000, Tuffs et al. 2004, Pierini et al. 2004), the far-infrared 
and sub-millimetre wavelengths probe its cool dust content, which 
is itself crucial to our understanding of star formation, since ap- 
proximately half of the energy ever radiated by stars has been ab- 
sorbed by dust and re-radiated at these wavelengths (e.g. Puget et 
al., 1996, Fixsen et al., 1998). The sub-millimetre region has been 
a difficult part of the electromagnetic spectrum in which to conduct 

* Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA 
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large galaxy surveys (e.g. Smail et al. 1997, Hughes et al. 1998, 
Eales et al. 1999). Previous ground-based sub-millimetre surveys 
had to be either pointed at pre-selected targets, or limited to rel- 
atively small regions of sky covering areas of < 1 deg 2 (Coppin 
et al., 2006, Weiss et al., 2009). The combined effects of the large 
negative fc-correction at these wavelengths, sensitivity and the steep 
number counts, have meant that the average 850 /im selected sub- 
millimetre galaxy is extremely luminous (10 12 — 10 13 Lq) and 
resides at high redshift (z ~ 2, e.g. Chapman et al. 2005). Few 
relatively local galaxies have been found in blind sub-mm surveys, 
due to the small local volumes probed in these surveys coupled with 
the observing wavelength targetting the faint Rayleigh- Jeans tail of 
the dust SED at low redshift. Our undertstanding of the local Uni- 
verse at sub-mm wavelengths has come so far from targetted sur- 
veys such as the SCUBA Local Universe Galaxy Survey (SLUGS, 
Dunne et al. 2000), which observed a sample of 184 IRAS- and 
optically-selected galaxies (Vlahakis, Dunne & Eales, 2005). Pre- 
selected galaxies in this way can lead to biases if there are classes 
of sub-mm emitting galaxies which are not bright at either optical 
or 60 mum wavelengths. The SLUGS surveys were also limited 
to very nearby galaxies and so could not address the question of 
evolution of sub-mm properties in the relatively recent past. 

With the advent of the PACS (Poglitsch et al. 2010) and SPIRE 
(Griffin et al. 2010) instruments aboard the ESA Herschel Space 
Observatory (Pilbratt et al. 2010), we now have our first oppor- 
tunity to survey a large area of sky at sub-mm wavelengths. The 
angular resolution and sensitivity of Herschel allow us to robustly 
determine the counterparts to thousands of local sub-millimetre se- 
lected galaxies across the whole electromagnetic spectrum, thus 
gaining invaluable insight into their physical processes. This pa- 
per uses a sample from the Herschel Astrophysical TeraHertz Large 
Area Survey (//-ATLAS: Eales et al. 2010) and presents fits to their 
UV-sub-mm SEDs. This is the first relatively local (z < 0.5) sub- 
mm selected sample for which such complete SED modelling has 
been performed. This work is based on only 3 percent of the final 
data-set, but is still large enough to provide a statistical study of the 
optical and IR properties of > 1000 250 ^m selected galaxies, and 
templates for SEDs which can be applied more widely. 
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Studies of the multi-wavelength properties of the relatively 
small number of galaxies detected in sub-millimetre surveys have 
been extensive (e.g. Swinbank et al. 2009). At high redshifts, 
galaxy star formation rates have been frequently estimated based 
on a single sub-millimetre flux measurement (e.g. at 850 /im), and 
a local template SED belonging to e.g. M82 or Arp 220 (e.g. Silva 
et al., 1998), chosen not because they are known to be representa- 
tive of the average sub-millimetre galaxy, but rather because they 
are comparatively well-studied. 

Another commonly-used method of describing far-infrared 
galaxy SEDs is to assume one or more components with modi- 
fied black-body (the so-called "grey-body") profiles. In these sim- 
ple parametrisations, the observed flux densities depend only on the 
temperatures (T) and dust emissivity index (/3), which may be ei- 
ther assumed or derived, depending on the available observations. 
Such simple grey-body profiles have been widely shown to broadly 
reproduce the sparsely-sampled far-IR SEDs of galaxies at all red- 
shifts (e.g. Dunne et al. 2000, Blain et al. 2002, Blain, Barnard & 
Chapman 2003, Pope et al. 2006, Kovacs et al. 2006, Dye et al., 
2010), although when the SED is sampled from 60/im to the sub- 
mm additional greybody components may be required to reproduce 
the observations (e.g. Dunne & Eales, 2001, Galametz et al., 201 1, 
Smith et al. 2010, Dale et al., 2012). 

Empirical templates have been created for use with sparsely- 
sampled far-infrared data based on observations of small samples of 
local galaxies with good coverage from mid- to far-IR wavelengths 
(e.g. Chary & Elbaz, 2001, Dale & Helou, 2002, Rieke et al. 2009). 
Selecting galaxies at shorter FIR wavelengths tends to favour those 
with substantial warm dust components, which may not be repre- 
sentative of populations selected at longer wavelengths with Her- 
schel and ground based sub-mm instruments. Several studies have 
found that sub-mm selected galaxies (so far mostly at higher red- 
shifts) may have colder dust than their local equivalents at similar 
far infrared luminosities (e.g. Pope et al. 2006, Coppin et al. 2008, 
Hwang et al., 2010, Magnelli et al., 2012). 

In this paper, we use a model that relies on energy balance - 
the idea that the energy absorbed by dust at ultra-violet and optical 
wavelengths must be re-radiated in the far-infrared - combined with 
a statistical fitting approach, to consistently model each galaxy's 
full SED, and gain robust constraints on the star formation activ- 
ity, stellar and dust content of 250 /im selected galaxies from H- 
ATLAS. 

In Section [2] we discuss the Herschel- ATLAS survey and the 
multi-wavelength data used in generating the catalogue, while in 
Section [3] we discuss the SED-fitting method used in the analy- 
ses which we present in Section [4] In Section [5] we compare our 
median SEDs with other templates available, and in Section [6] we 
present some conclusions based on our results for the population of 
sub-millimetre galaxies in general. Throughout this paper, we use 
a standard cosmology with Ho = 71kms~ 1 Mpc ~ , Qm = 0.27 
and 51a = 0.73, and an Initial Mass Function (IMF) from Chabrier 
(2003). 



2 CATALOGUE CONSTRUCTION 

The Herschel ATLAS (Eales et al. 2010) is the widest area sur- 
vey being conducted with Herschel covering 570 square degrees 
of sky in five FIR-sub-mm bands from 100-500 /im. One primary 
aim of //-ATLAS is to provide a census of dust and obscured star 
formation in the local Universe, with galaxies selected on the ba- 
sis of their dust mass for the first time. Our current study is based 



on the Herschel-KTLAS Science Demonstration Phase data (SDP) 
covering ~ 14 deg 2 centred on the 9hr GAMA field (Driver et 
al. 2011). The SPIRE and PACS map-making procedures are de- 
scribed in Pascale et al. (2011) and Ibar et al. (2010). From these 
maps a catalogue of sources which are ^ 5<r in any of the three 
SPIRE bands was produced using the MADX algorithm (Maddox 
et al. in prep) and described in detail in Rigby et al., (201 1). PACS 
sources were added to the catalogue based on the flux in apertures 
placed at the locations of SPIRE 250 /im sources. The catalogue we 
used for this sample is 250 /im selected and contains 6621 sources 
at > 5cr (though not all of these are detected in all other bands - 
see below and Table [TJ. The 5<r point source flux limits are 132, 
126, 32, 36 & 45 mjy in the 100 /im, 160 /im, 250 /im, 350 /im & 
500 /im bands, respectively (and including confusion), with beam 
sizes of approximately 9, 13, 18, 25 and 35 arcsec FWHM in the 
same five bands. 

A likelihood-ratio analysis (LR - Sutherland & Saunders, 
1992, Ciliegi et al., 2005) was performed to identify robust opti- 
cal counterparts to the sub-millimetre selected sources, using the 
SPIRE 250 /im channel and SDSS r-band positions down to a lim- 
iting magnitude of SDSS r modelmag = 22.4. The method used 
is described in detail in Smith et al. (2011), but to summarize; the 
LR method uses both positional and photometric information of 
both individual sources and of the population in general to quantify 
the reliability, R, (or equivalently, the probability) of an associa- 
tion between two sources. Star-galaxy separation was performed 
following a method similar to that in Baldry et al. (2010), and 
the LR calculations were applied to each population separately, as 
the 250 /im properties of stars and galaxies are quite different (see 
Smith et al., 201 1). For this study we have chosen a reliability limit 
of R ^ 0.8 which gives 2417 250 /im sources with reliable galaxy 
counterparts, and a contamination rate of less than five percent. We 
have also removed the five gravitationally lensed SMGs identified 
in Negrello et al. (2010) from the subsequent analysis. 

These data were combined with the Galaxy And Mass As- 
sembly (GAMA - Driver et al. 201 1) catalogue over the same field 
(Hill et al. 2011), which contains thousands of spectroscopic red- 
shifts, in addition to r-band defined aperture-matched photometry 
for 1402 of the reliable galaxy counterparts. Since we require well- 
sampled multi-wavelength SEDs to constrain the physical proper- 
ties of these 250 /im-selected galaxies, we base our analysis on this 
sub-sample of 1402 galaxies, with optical/near-infrared aperture- 
matched photometry. The r-defined aperture-matched photometry 
is based on pixel- and seeing-matched images derived from the 
Sloan Digital Sky Survey (SDSS - York et al., 2000) and the UK 
Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS - 
see e.g. Lawrence et al. 2007) in the ugrizY J H K bands. Addi- 
tional photometry in the Galaxy Evolution Explorer (GALEX) far- 
and near-ultraviolet channels has been included from the GALEX- 
GAMA survey (Seibert et al., in prep.), as have spectroscopic red- 
shifts from the GAMA, SDSS and 6dFGS surveys. There are 1052 
sources with spectroscopic redshifts and r-defined aperture pho- 
tometry. For the 350 sources with no spectroscopic redshifts avail- 
able, we adopt ANNz-derived (Collister & Lahav, 2004) neural 
network photometric redshifts from Smith et al. (201 1). 

All available photometry has been brought on to the AB mag- 
nitude system. Simulations have shown that the photometric errors 
estimated by S EXTRACTOR (Bertin & Arnouts, 1996) are underes- 
timated by a factor of four in the GAMA resampled images (Hill et 
al. 201 1); this factor was applied to the catalogue values. A further 
0.1 magnitude error was added in quadrature to all optical and NIR 
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Table 1. Table showing the coverage of the 250 ftm sources in our catalogue, detailing which sources are detected in which far-infrared bands, and detailing 
the number of sources with photometry from GALEX. The sensitivity limits indicate the properties of the input catalogues, note that in order for a source to 
have aperture photometry, we require a 5cr 250 ftm detection with an R ^ 0.8 counterpart from the upper row. 



Instrument: Herschel-SPTSE Herschel-PACS IRAS GALEX 

Band: 250 ftm 350 ftm 500 ftm 160 ftm 100 fim Both 60 /jm FUV NUV Both 



Catalogue Sensitivity: 


> 5cr 


> 


> 


> 5ct 


^ 5ct 


^ 5ct(x2) 


> 5cr 


> 5cr 


^ 5(7 


> 5ct(x2) 


N(detections) 


6621 


5346 


1717 


304 


151 


117 


34 








N(Galaxies, R > 0.8) 


2417 


1636 


344 


245 


142 


111 


34 








Aperture Photometry 


1402 


902 


170 


197 


116 


93 


24 


529 


726 


522 


^■spcc 


1095 


710 


139 


199 


1 18 


98 


33 








Aperture Photometry & z spec 


1052 


682 


128 


186 


108 


89 


24 


520 


700 


513 



photometry to account for the global uncertainties in the total flux 
measurements and calibrations between the various surveys. 

Additional photometry was compiled for known detections 
from the InfraRed Astronomical Satellite ( IRAS) Faint Source Cata- 
logue (Moshir et al., 1992, Wang & Rowan-Robinson, 2009), using 
a two arcsecond maximum match radius between cross-identified 
counterpart positions, including those sources with updated optical 
counterparts based on the higher-quality SPIRE and PACS images 
discussed in Smith et al. (2011). As noted in Table Q] there are 
34 sources with detections in the IRAS 60 /im band. The errors on 
the IRAS fluxes were assumed to be 20 per cent (including cali- 
bration error) and we included upper limits in the IRAS 12, 25 and 
60ftm channels based on the 5<j sensitivity limits given in Wang & 
Rowan-Robinson (2009). 

To reflect the uncertainty in the SPIRE and PACS photometric 
calibration, the SPIRE errors had a factor of 15 per cent added in 
quadrature to the catalogue values, and the PACS errors had 10 per 
cent and 20 per cent added in quadrature to the errors in the 100 and 
160 fim bands, respectively (e.g. Poglitsch et al. 2010, Griffin et al. 
2010). Although we only require 250 fim fluxes for the far-infrared 
selection, we include SPIRE fluxes in our input catalogue for each 
SPIRE band, irrespective of the signal to noise ratio in the 350 and 
500 ftm bands, provided that their measured flux is positive. For 
PACS, we only include those sources detected at ^ 5a significance 
in each band. This is due to residual 1 // noise in the current version 
of the PACS maps which limits the level at which we can extract 
reliable photometry at this current time. We plan to alleviate these 
problems in future releases. 

The number of sources with coverage in each of the far- 
infrared photometric bands, as well as the number of sources with 
aperture-matched photometry and spectroscopic redshifts are listed 
in Table Q] To summarize, our analysis will focus on a sub-sample 
of 1402 (i.e. 21 per cent) of the 5a 250 fim detections for which 
well-matched multi-wavelength data are available. 



2.1 The impact of selection effects 

Since the sample presented in this paper is not selected purely at 
250 fim, and relies on the identification of an optical counterpart 
brighter than 22.4 mag in the SDSS r-band data, it is necessary to 
consider the impact that this additional selection criterion might 
have on the results of this study. 

In the top panel of Figure Q] we show the variation in r-band 
magnitude for a representative template spectral energy distribu- 
tion based on the galaxies in our sample (see Section H31 for more 
details of how this SED template was derived). At each redshift be- 



ing considered, we fix the 250 fim flux of the template to the 5 a 
limit of our survey data (which corresponds to ~ 12.5 mag in the 
AB system), and convolve the template with the SDSS r-band fil- 
ter curve to determine the expected r-band magnitude that would be 
observed. We are then able to compare these values to the limits for 
the cross-identification (r mo d c i = 22.4 mag) and for the GAMA 
aperture-matched photometry (r po t ro = 20.5 mag), both of which 
are shown as the black labelled points and coloured lines in the 
upper panel. In the bottom panel of Figure [T] we compare the pre- 
dicted observed r-band magnitude as a function of redshift (solid 
line, with the surrounding shaded region corresponding to the un- 
certainty in the //-ATLAS template SED). The GAMA aperture- 
matched photometry limit is shown as a horizontal blue line. This 
shows that typical 250 ftm-selected galaxies should still be more 
than a magnitude brighter than the aperture-matched photometry 
limit even out at z = 0.5. 



We may still lose galaxies of certain types from the sam- 
ple at higher redshifts due to the optical limit, such as those with 
higher dust obscurations (i.e. higher 250 fim to r-band flux ratios). 
To investigate this further, in Figure [2] we show histograms of the 
(r— 250 fim) colour (in AB magnitudes) of those galaxies with re- 
liable counterparts in our sample, for different redshifts spanning 
0.0 < z < 0.5. Counterparts with aperture-matched photome- 
try in the GAMA catalogue (i.e. those galaxies to which we apply 
our SED fitting method) are shown as the filled grey histograms, 
while galaxies without are hatched red. The orange vertical dashed 
line indicates the (r— 250 fim) colour of sources at the sensitivity 
limits in both the r band and 250 ftm data, while the blue vertical 
dotted line shows the same colour for sources at the 250 fim flux 
limit and the nominal completeness limit for the GAMA photom- 
etry (20.5 mag in Petrosian magnitudes). It is clear from Figure [2] 
that our r-band selection criterion does not prevent us from obtain- 
ing a representative sample of 250 fim sources with a full spread of 
(r— 250 fim) colours, at least out to z < 0.35. At redshifts higher 
than z ~ 0.35 however, the r-band selection criterion does suggest 
that we are biased towards the lower obscuration sources (see also 
Dunne etal., 2011). 
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Figure 1. Top: expected magnitude of a median 250 /ira selected galaxy 
template (discussed in Section 14.31 normalised to the //-ATLAS 250 fim 
flux limit at redshifts between 0.0 < z < 0.5. The r-band limits for 
the optical cross-identification from Smith et al. (2011) and the GAMA 
aperture-matched photometry (Hill et al., 2011) are shown by the labelled 
blue and orange lines on the left hand side. The response curves of SPIRE 
at 250 fim and the SDSS r-band data are shown in the top corners. Bottom: 
the variation of the expected brightness of a median 250 (im selected galaxy 
template in the SDSS r band as a function of redshift. The aperture-matched 
photometry limit of r po t ro = 20.5 is shown by the blue horizontal line (i.e. 
the same colour as in the upper panel); the median 250 fim selected galaxy 
is at least a magnitude brighter than the r-band aperture-matched photome- 
try limit even out to z ~ 0.5. The shaded area corresponds to the range of 
SEDs in the //-ATLAS template SED. 



3 METHOD 

3.1 Spectral energy distribution modelling 

We use the model of da Cunha, Chariot & Elbaz (2008 - here- 
after DCE0£1) to interpret the panchromatic SEDs of the galaxies 
in our Herschel-ATLAS/IRAS/GAMA/GALEX data, set in terms of 
physical properties related to their star formation activity and dust 
content. This physically-motivated model relies on an energy bal- 
ance technique to interpret the (attenuated) stellar emission at ultra- 
violet, optical and near-infrared wavelengths consistently with the 
dust emission at mid/far-infrared and sub-millimetre wavelengths. 
Therefore this model is ideal to interpret the multi-wavelength ob- 
servations available for this sample of galaxies. Here we briefly 
summarize the main features of this model; for more details we 
refer to the exhaustive description provided in DCE08. 

The dust-free ultraviolet to near-infrared emission from stel- 
lar populations in galaxies is computed using the latest version of 
the Bruzual & Chariot (2003) stellar population synthesis models 
(Bruzual 2007; Chariot & Bruzual, in prep.). The attenuation of 
starlight by dust is described by the two-component prescription 
of Chariot & Fall (2000), which also provides the total energy ab- 
sorbed by dust in the birth clouds (i.e. molecular clouds where stars 
form) and in the ambient (i.e. diffuse) interstellar medium (ISM). 
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Figure 2. The (3 — 250 fim) colours (in AB magnitudes) of the galaxies in 
our sample, for different bins of redshift from 0.0 < 2 < 0.10 (top), to 
0.35 < z < 0.40 (bottom). The grey histograms represent those galaxies 
for which we have aperture-matched photometry from the u—K bands from 
the GAMA survey, while the dashed red histograms show the colours of 
those galaxies for which this aperture-matched photometry is not available 
(these are fainter than r po t r o = 20.5 or below 8 = —1.0° declination due 
to the GAMA magnitude limit, and to the peculiarities of the GAMA///- 
ATLAS survey overlap, respectively). The number of sources falling into 
each of these categories is shown in the top right hand corner of each panel. 
The vertical dashed line indicates the (r— 250 ^tm) colour corresponding 
to the sensitivity limits of our r-band and 250 fim catalogues, while the 
vertical dotted line indicates the colour defined by the magnitude limit for 
the GAMA aperture-matched photometry relative to the 250 fim flux limit. 
It is clear that we sample a representative range of 250 fim selected galaxy 
colours out to z fa 0.35, but beyond this we are only sensitive to the less 
obscured part of the population. 



The spectral distribution of the energy re-radiated by dust at infra- 
red and sub-millimetre wavelengths is then computed by assuming 
that the energy re-radiated by dust in the birth clouds and diffuse 
ISM is equal to the energy absorbed, and that starlight is the only 
significant source of heating (i.e. that there is no AGN contribu- 
tion). In stellar birth clouds, the dust emission is described as a sum 
of three components: polycyclic aromatic hydrocarbons (PAHs), 
hot mid-infrared continuum and warm dust in thermal equilibrium 
with temperature in the range 30 - 60 K. In the ambient ISM, the 
emission by dust is described using these three components (whose 
relative proportions are fixed for simplicity), plus a component of 
cold dust in thermal equilibrium with temperature in the range 15 
- 25 K. The prior distribution of both the warm and cold dust tem- 
peratures is flat, such that all temperatures between the bounds of 
the prior have the same probability^ The hot and cold dust com- 
ponents of the spectrum are assumed to be optically thin, and are 
described in the same way as in DCE08 using modified greybody 
template spectra, oc k\B\(T), with emissivity index j3 — 1.5 and 
2.0 for the warm and cold components, respectively, and dust mass 
absorption coefficient approximated as a power law, such that 



The da Cunha, Chariot & Elbaz (2008) models are publicly available as 
a user-friendly model package MAGPHYS at: www.iap.fr/magphys. 



2 Our choice of prior distribution for the cold dust temperature is discussed 
in more detail in Appendix lAl 
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KA CC A P , (1) 

with the normalisation defined such that KssOfim = 0.077 kg -1 m 2 
as in Dunne et al., 2000. 

The simplicity and versatility of the DCE08 model make it 
ideal to interpret our rich multi-wavelength data set, as it allows us 
to derive statistical constraints (including probability density func- 
tions, hereafter PDFs) for several physical properties of the galaxies 
(such as star formation rate, stellar mass, dust attenuation, dust lu- 
minosity measured between 3-1000 ^im, dust temperature and dust 
mass), from the consistent modelling of their observed ultraviolet 
to sub-millimetre spectral energy distributions. To do so, we adopt 
the Bayesian approach used in DCE08 (see also da Cunha et al. 
2010, hereafter dCIO). 

We use two stochastic libraries of models, as described in 
DCE08; the first contains 25,000 stellar population models, in- 
cluding a wide range of star formation histories, metallicities and 
dust attenuations, while the second consists of 50,000 dust emis- 
sion models including a large range of dust temperatures and frac- 
tional contributions of PAHs, hot mid-infrared continuum, warm 
dust and cold dust to the total infrared luminosity. These two li- 
braries are combined by associating models with similar values 
of / M = ijf,^t /^dust ( me fraction of total dust luminosity con- 
tributed by the diffuse ISM), which are scaled to the same total dust 
luminosity LdSst • F° r eacn combined model spectrum, we compute 
a library of synthetic photometry making use of the filter transmis- 
sion curves for the same photometric bands as our observations, at 
intervals of dz = 0.01. We do not consider redshift to be a free pa- 
rameter in our SED fits, and use either the spectroscopic redshift or 
the best-fit photometric redshift from Smith et al. (2011) for these 
purposes. For ease of reference, some of the output parameters to 
which we refer in this analysis are summarized in Table|2] 

3.2 Spectral fits 

We use the MAGPHYS model to fit the observed SEDs for the 1402 
Herschel-ATLAS sources with robust (i.e. Reliability > 0.8) coun- 
terparts and matched-aperture photometry in the GAMA catalogue. 
For each galaxy we compared the observed ultraviolet to sub-mm 
fluxes to the predicted fluxes of every pair of models in the stochas- 
tic libraries satisfying the energy balance criterion, by computing 
the x goodness-of-fit parameter for each model. This allowed us 
to build the PDF of any given physical parameter for the observed 
galaxy by weighting the value of that parameter in each model by 
the probability exp(— x 2 /2). We then determine the median value 
of each PDF, corresponding to our best estimate for each param- 
eter. We also determine an associated uncertainty, which corre- 
sponds to the 16th & 84th percentiles of the PDF. We note that the 
PDFs generated in this analysis are marginalised; this is particularly 
important since by definition they include parameter uncertainties 
due to e.g. the co-variances between parameters in the model. In 
what follows, the values of the physical properties of the galaxies 
mentioned refer to the median values of the PDF unless explicitly 
stated otherwise. We also create 'stacked PDFs' when discussing 
the parameter properties of samples of sources, this procedure is 
described fully in Appendix Al but is designed to give our best 
estimate of the distribution of parameter values for sources in a 
sample convolved with our ability to constrain them. In addition to 
the PDFs for each model parameter, we also obtained the best-fit 
model SED for each galaxy, which is the model that minimizes \ 2 . 
It is important to determine whether or not the best-fit model 



reasonably reproduces our observed data, which are not uniform 
across the sample. Since neighbouring photometric bands are not 
independent of one another, we conducted several sets of simula- 
tions, designed to empirically estimate the variation of the number 
of degrees of freedom in our spectral fits, as a function of the num- 
ber of bands of input photometry available for a particular galaxy. 
The details of these simulations are presented in detail in Appendix 
|B1 To summarize, we determined a 99 percent confidence interval 
on x 2 , which depends solely on the number of photometric detec- 
tions for a particular source, such that if the derived value of \ 2 is 
higher than the upper bound of the interval, there is only a < 1 per 
cent chance that the galaxy is well described by our model, and is 
then removed from our sample. 

We find that ~ 92 per cent of our sample are well-described 
by our model. The galaxies with larger \ 2 values than our limit may 
have problems with photometry, contributions from AGN compo- 
nents, be lensed systems, or have catastrophic photometric redshift 
errors etc. Inspection of the bad fits (~8 per cent) reveals that 
the vast majority are due to serious problems with the aperture- 
matched photometry (e.g. catastrophic failures in multiple bands), 
while at least one is a QSO, and there are two possible lensed ob- 
jects with far-IR colours not consistent with the redshifts of their 
R ^ 0.80 counterparts, similar to those discussed in Negrello et al. 
(2010). 

We note that there are 320 galaxies in our sample which rely 
on photometric redshift estimates, of which 126 lie at z < 0.35, 
where our sample is thought to be representative of the broader 
250 j-im selected population. Nine per cent of those galaxies rely- 
ing on photometric redshifts have \ 2 values outside the range of 
acceptable values, as compared with seven per cent of those galax- 
ies with spectroscopic redshifts, suggesting that sources with catas- 
trophic photometric redshift errors do not make up the majority of 
unreliable SED fits. 

Finally, we note that dCIO conducted a series of tests of this 
SED-fitting model using a local galaxy sample detected with SDSS 
and IRAS, and demonstrated the general robustness of the technique 
to the effects of inclination. They showed that while weak inclina- 
tion effects may be present (using the ratio of the apparent major 
and minor axes of each particular galaxy as a proxy for inclina- 
tion), they do not dominate the dispersion in estimates of galaxy 
properties. 



4 RESULTS 

Here we present the results of our SED-fitting analysis. For each 
galaxy we not only determine best-fit SEDs (Figure[3j, but also the 
PDFs for each parameter (Figure[4]l. The parameters that we focus 
on are those shown in Figure|4](namely M dust , L d ust, M sta rs, f^, 
SFR, sSFR, Mdust/Mstars and M du at/£dust), since we are inter- 
ested in investigating the star formation activity and dust mass of 
normal, star-forming galaxies. Our sample represents almost an or- 
der of magnitude's improvement upon the largest sub-mm selected 
samples until now, even though the H- ATLAS SDP data comprise 
only ~3 per cent of the total eventual //-ATLAS data set. 

4.1 Importance of the available infrared observations 

Given the inhomogeneous set of data which we have compiled for 
the //-ATLAS sample, it is important to have some understanding 
of the sensitivity of the derived physical parameters to the absence/ 
presence of data at certain wavelengths. Some tests of this nature 
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Table 2. Summary of the SED model parameters used in this paper. For a more detailed description of each parameter, see da Cunha, Chariot & Elbaz (2008). 



Parameter Definition 



u 

sSFR 
SFR 

Mdust 
Ldust 
A/stars 
jfcold 



Fraction of total dust luminosity contributed by the diffuse interstellar medium 

Specific star formation rate is defined as the star formation rate per unit stellar mass, averaged over the last 0. 1 Gyr, units of yr 
Star formation rate averaged over the last 0.1 Gyr in Mq yr -1 
Dust mass in solar units, Mq 

Dust luminosity integrated between 3 & 1000 fim, in units of Lq 
Galaxy stellar mass in units of Mq 

Temperature of the cold dust in thermal equilibrium in the diffuse ISM 



e io 10 




0.1 



1.0 10.0 

Rest-frame Wavelength (|j.m) 



100.0 



1000.0 



Figure 3. An example best fit SED for one of the galaxies in our sample (HATLAS J090713. 1-000322). This SED is derived based on our compiled data 
from SPIRE and PACS, plus upper limits at 12, 25 & 60 from IRAS, and aperture-matched photometry from u-K bands as well as photometry in the 
GALEX FUV and NUV bands. The data points are shown as the black crosses, with the error bars as discussed in Section[2] The upper limits are displayed 
as down-pointing arrows. The best fit total SED is plotted in grey (thick), with the stellar component in red and the corresponding dust model in blue. The 
unattemuated (i.e. dust-free) model stellar SED is shown in green. The most uncertain area of the SED is the mid-infrared, which is constrained only weakly 
by the IRAS upper limits, with additional constraint coming from the energy balance criterion. The energy balance technique of da Cunha et al. (2010) enables 
us to only combine optical/near infrared SEDs that are physically consistent with the sub-millimetre SED. 



were performed by DCE08 and dCIO, and demonstrated that the 
results derived were generally robust, although certain parameters 
(e.g. sSFR) are better constrained when UV data are included, for 
example. Our concerns for this study relate to determining the re- 
liability of the estimates of dust luminosity, SFR and mass trends. 
Given that not all galaxies in the sample have data spanning the 
peak of their rest-frame far-infrared SED (from IRAS or PACS) we 
need to assess the impact of this heterogeneity on our results. We 
have performed three tests. The first test was to determine the re- 
liability of our method when only IRAS data between 60-100 pm 
are used to constrain the far-infrared SED; secondly, we attempted 
to determine the influence of incomplete PACS data on our results, 
and finally, we attempted to probe the reliability of our results in the 
absence of mid-IR data, which are not available over the //-ATLAS 
SDP field at the time of writing. 



4.1.1 Comparison with IRAS-selected samples 



In Section POl we compare the star formation activity and dust con- 
tent of //-ATLAS 250 /im-selected galaxies with those of a previ- 
ous sample of local, star-forming galaxies selected at 60 /im with 
IRAS (da Cunha et al., 2010). To compare these two samples, we 
need to understand possible differences/biases in the derived phys- 



ical parameters that may arise from the different selection of the 
samples. Therefore, in this section, we investigate the effects of in- 
cluding SPIRE data in the SED fitting for a sub-sample 250 pm 
selected galaxies for which we also have IRAS data. This allows us 
to assess whether any difference between our results and those of 
dCIO are due solely to the lack of SPIRE data for that study, or if 
they are rather due to effects of selection. 

We take a sub-sample of the //-ATLAS 250 /im-selected 
galaxies which are also detected by IRAS at 60 /im, and applied our 
fitting procedure twice; once including all available data, and a sec- 
ond time omitting all data longward of the PACS 100 pm band. The 
PACS 100 ^m data were included for both sets of fitting to ensure 
that our results are as comparable as possible with dCIO, since their 
galaxy sample was bright enough to be well-detected with IRAS 
at 100 pm. There are 18 //M5-detected galaxies with aperture- 
matched u- to A'-band photometry which are well-described by our 
models. For these galaxies fitted both with and without the 160- 
500 pm data, we compare the medians of the stacked PDFs for a 
given parameter in each run. We estimate the uncertainty associated 
with each bin in a stacked PDF according to the 16th and 84th per- 
centiles of the cumulative frequency distribution of values in each 
PDF bin (this method is discussed in greater detail in Appendix lAt. 
These stacked PDF comparisons and their associated uncertainties 
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Figure 4. Probability density functions (PDFs) for the galaxy shown in figure [3] (HATLAS J090713. 1-000322); the dust mass, dust luminosity, stellar mass, 
/^ (the fraction of the total dust luminosity contributed by the diffuse ISM), the [specific] star formation rate, dust to stellar mass ratio, and the ratio of dust 
mass to dust luminosity. Our results suggest that this galaxy is a massive luminous infrared galaxy (LIRG). Using these PDFs we can estimate not only the 
median values for each parameter, but also the 16th & 84th percentiles, corresponding to the confidence range of our estimates. 
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Figure 5. Stacked probability density functions and errors for sources in our sample detected by IRAS at 60 /im, derived using the full IRAS, PACS, & SPIRE 
data set (red histograms), and using only the IRAS 60 (im and PACS 100 fim data (blue histograms). The median values of each stacked PDF are indicated by 
the vertical dashed lines. The increased uncertainty in dust mass when the A > 100 fim data are not included in the fits is apparent in the wider range of values 
in the blue (i.e. no Herschel) M dust PDF. Stellar masses, / M , dust luminosity and star formation rates are much less strongly influenced by the inclusion of 
longer wavelength data i.e. they are well constrained without data beyond 100 fjm. 



are displayed in Figure[5] with the full data set PDFs shown in red, 
and the IRAS-onXy results shown in blue. 

In general, removing the A ^ 160 /im data from the SED fit- 
ting for those sources detected by IRAS results in small variations 
of the median values of the stacked PDFs for each population (dot- 
ted vertical lines in Figure[5}, but the changes are always less than 
the confidence interval derived according to the 16th and 84th per- 
centiles of the stacked PDFs. The most noticeable effect of remov- 
ing the A ^ 160 data is the considerably larger uncertainty in 



Mduat; this is not surprising since the cold dust component domi- 
nates the total dust mass, and constraints on this component come 
primarily from the longer wavelength SPIRE data. There is a ten- 
dency for the model to overestimate the dust mass when using only 
IRAS data (see the tail to high Md us t values in Figure|5]l, as the lack 
of SPIRE data allows the model to add in more cold dust without 
any strong constraint from the energy balance (since the contribu- 
tion to the overall Ld us t varies oc T 4 + 8 5 warmer dust far outweighs 
colder dust in its effect on .Ldust 

). The quantities M duat /L dust , 
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Figure 6. S250 colour as a function of 250 fim flux density, binned 
in redshift according to the colours as shown in the legend. At the fainter 
250 (im fluxes our PACS data detect only the warmer galaxies, as shown by 
the large number of lower limits on the S250 / Sieo colour at these fluxes. 
At S"25() > 120 mJy we sample the full range of colours (and therefore 
temperatures) in our 250 fim selected sample, so it is these galaxies that we 
use to study the impact of missing PACS data for the majority of our sources 
on the derived parameters. The vertical dashed line shows S250 = 120 mJy, 
the 250 fim flux limit for the PACS-complete sub-sample. 

and Mdust/Afstars are much better constrained when including the 
longer wavelength data, for the same reasons discussed for the dust 
mass estimates above. The other parameters, Ldust, Mstars, SFR, 
sSFR and / M are all comparable in both samples. 

4.1.2 The effect of lacking PACS observations 

Whilst there are 1289 5a SPIRE sources with reliable optical coun- 
terparts in the //-ATLAS survey and good SED fits, we have 5<r 
PACS detections for only 207 and must rely on upper limits for the 
remaining sources. It is important to understand the effect of miss- 
ing PACS data on our results, and so we now investigate the impact 
on our SED fits when we omit the PACS data for a sub-sample 
of PACS-detected galaxies drawn from our main sample. To deter- 
mine a representative sample, we show in Figure[6]the S'250/5'160 
colour as a function of 250 fim flux, binned in redshift. The PACS 
data are not deep enough to probe colder colours at fainter fluxes 
or higher redshifts, but do sample the full range of colours above 
S250 ~ 120 mly, at least for z < 0.2. To assess the impact of 
missing PACS data on the results of our SED fitting, we use these 
59 sources with 160/im PACS detections and 52so~ 120 mly at 
z < 0.2, to check how the SED parameters change when the PACS 
data are removed from the fitting^. We will refer to this sample as 
the "PACS-complete" sample. 

In Figure[7] we show a comparison of the stacked PDFs when 
the PACS data are included and excluded from the fitting for the 
"PACS-complete" sample. This is not quite the same effect which 
will apply to the full sample, as there we do have some PACS in- 
formation (in the form of upper limits), while in this test we have 
no information when the PACS data are removed. This makes this 
comparison conservative, in the sense that the effects on the full 
sample are going to be no larger than the worst-case scenario we 
study here. 

Removing the PACS information from the SED fitting causes 
the average estimate of the specific star formation rate to decrease 

3 the equivalent sample for the PACS 100 fim channel, which requires 
S2S0~ 200 mJy, contains no additional sources above those selected here 



by ~ 0.1 dex. Estimates of the stellar mass are, unsurprisingly, 
barely altered, while / M and Ldust are also robust in a sample av- 
erage. Dust masses of the population are overestimated by approx- 
imately ~ 0.06 dex. The median offset for each parameter in the 
absence of PACS data is tabulated in Table[3] 

Looking at the changes in globally averaged parameters in this 
way is reassuring, however, it is important to check that this is not 
masking a potential correlation of a bias in one parameter as a func- 
tion of another. For example, we may overestimate Mtust at low 
values of Ldust and underestimate at high values - i.e. our esti- 
mates may be skewed - but still have an average offset consistent 
with zero bias. As we will next investigate trends of one parame- 
ter against others, and later bin SEDs by parameter, we must con- 
sider these effects now. We focus on these issues in detail in Ap- 
pendix|C] though to summarise, we find that the lack of PACS data 
does not introduce bias in any parameters as a function of redshift, 
M a tars, or f,j.. The same is generally true of our estimated Ldust 
and Mdust, though these parameters possibly show weak correla- 
tions (i.e. they may be skewed). In the absence of PACS data, Ldust 
may be under-estimated at large Mdust, and over-estimated toward 
lower dust masses, though the offset averaged over the whole popu- 
lation is small. Mdust shows weak bias with Mdust ACS ; m that our 
PACS-free estimates of A/dust are slightly high for large Mdust, 
though the correlation shows considerable scatter and the overall 
offset across all values is small. Our estimates of SFR and sSFR 
show larger scatter than the other parameters, reflected in the larger 
error bars on A(s)SFR, though it is difficult to discern any skewed 
bias in the derived values. We will discuss the impact of possible 
bias with these parameters later in this paper. 

It is worth noting that because of the good multi-wavelength 
coverage and lack of temperature/colour bias in the S250 > 
120 mJy PACS-complete sample (Figure [6j, we can use it to de- 
termine our best estimate of the dust temperature in low-redshift 
250 /im-selected galaxies. We determine a median likelihood es- 
timate of T 0O id = 20.4 ± 0.4 K from the DCE08 model, and an 
isothermal value from simple grey-body fitting to the FIR data 
of Tg r ey = 26.1 ± 3.5 K (assuming /3 = 1.5 for compari- 
son with literature values). These estimates are colder than values 
in the literature pre-dating Herschel; for example, SLUGS found 
Tgrey = 36 ± 5 K, /3 = 1.3, for the //MS-selected sample (Dunne 
et al., 2000) and T grcy = 31.6 ± 0.6 K, for the optically-selected 
sample (Vlahakis, Dunne & Eales, 2005). Our temperatures are 
comparable to those observed in galaxy samples selected at longer 
wavelengths but including data from 24-160 using e.g. BLAST 
(26 ± 5 K, /3 = 1.5, Dye et al., 2009) or selected in the K-band and 
observed with Herschel (e.g. T groy ~ 20 K, f3 — 2.0, Boselli et al., 
2010, T grcy « 25.8 K for spirals in Skibba et al. 201 1, or w 23 K, 
P ~ 1.5, Dale et al., 2012). Multiple-component modified black- 
body SED fits have long noted the presence of substantial cold dust 
components with lower temperatures, consistent with our findings 
(e.g. Dunne & Eales 2001; Contursi et al., 2001; Vlahakis, Dunne 
& Eales 2005; Clements, Dunne & Eales, 2011, Rowan-Robinson 
et al., 2011; Galametz et al, 2011, Planck Collaboration, 2011). 

4.1.3 The effect of missing mid-IR observations 

At the time of writing, mid-infrared data over the //-ATLAS fields 
are unavailable (the first public data release of the Wide-Field Infra- 
red Survey Explorer - WISE, Wright et al., 2010 - survey does not 
include the //-ATLAS SDP field), we investigate the effects of not 
having such data on the PDFs by applying the same fitting proce- 
dure to a sample of galaxies selected at 250 in data from the 
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Figure 7. Stacked PDFs comparing the dust temperatures, dust mass and luminosity, stellar mass, , specific star formation rate, dust mass to luminosity 
ratio and dust to stellar mass ratio for the representative PACS-complete subsample. The results including all available data are shown by the red histograms, 
while those results derived neglecting the PACS data are shown in blue. The vertical dotted lines indicate the median values of the stacked PDFs. We find that 
fits lacking PACS data have SSFR underestimated by ~ O.f dex, and Af(j ust overestimated by ~ 0.06 dex. Estimates of the stellar mass and / M are largely 
unaffected by omitting the PACS data, while the effects on the M^ ust / L^^rnid M^ ua ^/M a ^ aIs e^e consistent wth the changes to their input parameters. 



Table 3. The effects of omitting PACS data from the SED fitting results for 
those sources in the "PACS-complete" sample. These offsets are also shown 
in Figure [7] Note that all values are consistent within the uncertainties de- 
rived from the 16 th and 84 th percentiles of the stacked PDFs. Furthermore, 
and as can be seen in figure |C1| of appendix [C] the offsets for individual 
galaxies are almost always consistent with zero once the errors are taken in 
to account. 



Parameter 


Offset 
(Best - no PACS) 


l°gl0 ^dust 


0.03 dex 


l°glO M dust 


-0.06 dex 


log 10 M stars 


0.00 dex 


u 


-0.01 


sSFR 


-0.10 dex 


SFR 


-0.10 dex 


A^dustAkdust 


-0.1 5 dex 


Mjust/Afstars 


-0.04 dex 



Balloon-bome Large Aperture Sub-millimetre Telescope (BLAST 
- Devlin et al., 2009). We derive multiwavelength SEDs and PDFs 
for a sample of 14 BLAST galaxies in the Extended Chandra Deep 
Field South (ECDFS, Lehmer et al., 2005) which have ^5 o de- 
tections in all three BLAST bands (250, 300, 500 /urn, Dye et al., 
2009), and spectroscopic redshifts from Eales et al. (2009), with 
additional photometry in the GALEX (FUV and NUV - Morrissey 
et al., 2007), optical {ugriz), 2MASS J and K s (Skrutskie et al., 
2006), MIR (3.6-8.0, 24, 70 /xm), as well as far-infrared (160 urn) 
bands from the Spitzer Space Telescope (Lonsdale et al., 2003). 

In Figure[8]we show the stacked probability density functions 
for the same parameters as in figure |7J derived for these BLAST 
galaxies. The red histograms show the PDFs determined when we 
include the complete data-set, while the blue histograms show the 
PDFs derived in the absence of the mid-infrared data. The similar- 
ity between these two sets of histograms, and the absence of bias 



between them, suggests that our estimates of the dust mass, lumi- 
nosity and SFR are robust to the absence of mid-infrared data in 
our wider //-ATLAS data set, although these tests have necessar- 
ily only been applied for a small number of sources. The mid-IR 
accounts for only a small fraction of the total infrared emission, 
and the similarity of the PDFs highlights the power of the energy 
balance criterion in constraining the dust luminosity even in the ab- 
sence of mid-infrared data. The detailed shape of the SED in the 
mid-IR is clearly not well defined for our sample and constraints 
on this can only come from comparison with mid-IR data, e.g. from 
WISE. 



4.2 The properties of Sub-millimetre selected galaxies in 
//-ATLAS 

4.2.1 First results, and comparison with previous studies 

By stacking the PDFs for galaxies well described by our model, we 
determine a median dust luminosity for our whole sample of 5.6 x 
1O 1O L0, placing the average //-ATLAS galaxy's luminosity just 
below what would traditionally have been considered a luminous 
infrared galaxy (LIRG). Figure [9] shows the dust luminosity as a 
function of redshift for //-ATLAS (red) and it is important to note 
that //-ATLAS traces typical star forming spirals (with log Ld us t < 
11.0 Lq) out to much higher redshifts (z ~ 0.35) than was possible 
with IRAS (z < 0.05). Our SED-fitting results indicate that the 
star formation rate of the average low-redshift //-ATLAS galaxy is 
~ 3.4 Mq yr" 1 , and that the median dust to stellar mass ratio is 
~ 0.4%. 

In FigureQJJ] we present histograms of the results of our SED 
fitting, and those from dCIO, who applied the same method to 
a GAL£X-SDSS-2MASS-/flAS data set and included photometry 
only up to 100 pm. The different selection criteria of the dCIO sam- 
ple produces a different redshift distribution from the //-ATLAS 
sample, however the stellar mass and dust luminosity distribu- 
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Figure 8. A comparison between the stacked PDFs derived by fitting 14 galaxies selected from the BLAST survey of the E-CDFS. In red we show histograms 
of the results obtained when the mid-infrared data (i.e. those at 3.6-8.0, 24 and 70 fim) are included in the fitting, whilst the blue histograms show the results in 
the absence of these data. The similarity of these two sets of PDFs indicates that our estimates of these parameters are not biased by the absence of mid-infrared 
coverage in our data set, and consequently that the absence of such data does not limit the validity of our results. Once more, the median values of the PDFs 
are indicated by the vertical dotted lines. 




Table 4. Comparison between the median properties of the dCIO IRAS- 
selected sample and the //-ATLAS 250 /im-selected sample. 
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Figure 9. The dust luminosity of Z/ersc/ieZ-ATLAS-selected galaxies with 
reliable counterparts, and reliable SED fits to our compiled photometry, as 
a function of redshift, is shown in red. The median redshift is z = 0.239 
and the median dust luminosity of our full sample is 5.6 X 10 10 Lq. We 
also include data points from dCIO for comparison, in black, highlighting 
the different redshift and luminosity properties of each sample. The hori- 
zontal dashed lines delineate those galaxies with luminosities in the LIRG 
and ULIRG categories, at 10 11 and 10 12 L0, respectively. We also show 
histograms detailing the redshift distribution of each sample (top), in which 
the relative heights of the two histograms are arbitrary to bring out the con- 
trasting properties of the two samples. On the right we show that the dust 
luminosity distributions for the two samples are comparable. 



tions are similar, and comparison of specific parameters can still 
be instructive. In Table [4] we provide results for the dCIO sample 
(z < 0.16), and for our sample of 250 fim selected galaxies (lim- 
ited to z < 0.35). 

On average, the galaxies selected at 250 /im have slightly 
higher specific star formation rates than those in dCIO, by 
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-2.381 
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0.57 
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0.1 dex. Though the dust luminosities and stellar 
masses are roughly the same, the dust masses in the two samples 
differ by 0.3 dex, with the //-ATLAS galaxies being more dusty. 
This translates into higher "specific" dust masses for //-ATLAS 
galaxies (higher M dust /Ldust and M dust /Af s tars)- Part of this dif- 
ference is due to the evolution in dust masses with redshift (Dunne 
et al. 201 1), as the //-ATLAS sample probes a higher redshift range. 
However, this is not the whole story as we will discuss in the next 
section. 

In figure[T0] we compare the stacked PDFs of several physical 
parameters for the //M5-selected sample of dCIO (blue) and the 
//-ATLAS 250 ^m-selected sample (red). The dCIO distribution of 
AfdustALdust values is considerably narrower than that based on 
the results of this study. This may be because Herschel's selection 
at 250 /im is intrinsically more sensitive to a range of M duat /L d ust 
values, as it is not solely sensitive to the warmer dust but also picks 
up the cold, dusty galaxies which have large M dust /I/ dus t but are 
not necessarily warm enough for IRAS to detect. 
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Figure 10. Comparisons of the stacked PDFs from //-ATLAS SED fitting (red histograms) with those of dCIO (blue histograms). The red dotted histograms 
indicate the distribution of the median likelihood values for //-ATLAS (i.e. the 50th percentiles of the individual PDFs, arbitrarily renonnalised for ease of 
comparison), which are narrower than the stacked PDFs since they do not reflect the uncertainty in the derived values. The histograms compare parameter 
distributions for the z < 0.35 sample which are well described by our model, with the results of dCIO. We note that the dCIO results do not include PDFs 
for SFR, M^ ust /L^ ust or Md us t/A^stars, and so we include the quotient of the median values of the relevant PDFs for the dCIO sample shown as thick blue 
dotted histograms. The stacked results suggest that galaxies selected in //-ATLAS contain on average a factor ~ 2 more dust than those galaxies in the dCIO 
sample, although the luminosities are broadly comparable. 



4.2.2 Star formation and dust in 250 fim-selected galaxies 

We now explore the star formation activity and dust properties of H- 
ATLAS galaxies. In Figure [TT] we plot relationships for three dif- 
ferent galaxy samples; the results of this study (left), and the dCIO/ 
DCE08 samples (right, orange and blue squares, respectively). 
Each individual galaxy is shown in grey, and we split our sample 
in bins of redshift, with colours as shown in the legend. The posi- 
tions of the error bars correspond to the mean of the galaxies in that 
redshift bin, while the size of the error bar represents the standard 
deviation of the derived values within that redshift bin. The best 
fit relation between SFR and dust mass from dCIO (dashed line) 
appears to trace the low dust-mass edge of the 250 /im selected 
population. The slope of the relationship for //-ATLAS appears to 
be the same, but the //-ATLAS sources are offset such that they 
have larger dust masses for a given SFR compared to IRAS selected 
galaxies. 

It is unlikely that this offset is a result of biases in the fit- 
ting given the difference far-infrared coverage for the two samples. 
As we showed in Section |4.1.1| for the IRAS sample, we do not 
expect A/dust to be under-estimated (if anything, the converse ap- 
plies) due to the lack of data at wavelengths longer than 100 /im. 
For //-ATLAS sources without PACS data, there was a small ten- 
dency to overestimate the dust mass at the highest values of SFR, 
but not at a level which could explain this offset which is present 
at all SFR values. This effect is also not due to redshift differences 
between the two samples, as the error bars show, increasing red- 
shift drives a given galaxy along the plotted slope and not away 
from it. Comparing only sources within the same redshift range 
(i.e. if we limit the //-ATLAS results to z < 0.16), we still find 
an increase in dust mass per unit star formation rate for the H- 
ATLAS sample, relative to the DCE08 and dCIO samples. At low 
values of dust mass, the //-ATLAS sample includes few sources 



with log 10 (Md us t/ M©) < 6.5. This is due to the flux limit in H- 
ATLAS combined with the small survey area in SDR The "missing 
//-ATLAS sources" following the IRAS and SINGS points at the 
lower left of the plot are simply below the detection threshold of 
the //-ATLAS SDP sample. 

The lack of IRAS sources with high A/d us t/SFR may reflect 
the lack of sensitivity of IRAS to cold dust; galaxies with low SFR 
and larger masses of dust would tend to have colder dust temper- 
atures and therefore be absent from IRAS selected samples. This 
comparison suggests that IRAS preferentially selects those galaxies 
with the highest SFR per unit mass of dust, since more star forma- 
tion for a given mass of dust will result in stronger dust heating and 
more emission at the IRAS selection wavelength of 60 /im. For the 
DCE08 SINGS sample (squares), which are not selected from a flux 
limited FIR survey, the trend is similar to that seen in //-ATLAS but 
extended to lower values of dust mass. 

In the bottom panel of Figure [TT] we also plot the ratio of the 
dust to stellar mass as a function of specific star formation rate. 
This relationship was first noted by dCIO, and their data are plotted 
in the right panel of Figure [TT] Compared to the sample in dCIO, 
the //-ATLAS galaxies appear to have higher specific dust content 
(relative to stellar mass) for a given specific star formation rate. 
The reason for this difference is likely to be the same as that in the 
upper relationship between Mdust and SFR since the stellar mass 
distributions of the two samples are very similar. 

We overlay the predictions of dust evolution tracks from the 
chemical evolution models of Calura et al. (2008), which are 
based on the model in Dwek (1998) and follow the build-up of 
heavy elements and dust formed in low-and-intermediate mass stars 
(LIMS) during their AGB phase and in both Type la and Type II 
supernovae. The upper black curve shows the evolution of dust 
from a dwarf-irregular galaxy with continuous star formation. The 
solar-neighbourhood model, which reproduces the properties of the 
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Figure 11. Comparisons between //-ATLAS 250^m-selected galaxies and other local galaxy surveys. Top Row: the variation of the total dust mass as a 
function of the star formation rate for the three samples of galaxies; on the left, we show the //-ATLAS 250^m-selected sample with at least least two far-IR 
detections in our data set are shown in grey. The coloured error bars indicate the mean positions of those galaxies in the redshift range shown by the colours, 
with size equal to the standard deviation of the values in each bin. On the right we present the same results for the GALEX-SDSS-2MASS-//14S sample from 
dCK) in orange circles, with the SINGS galaxy sample from dC08 in blue squares. Overplotted in each frame is the best fit relationship from da Cunha et al. 
2010, which is offset with respect to the results of our best fit SEDs. Bottom row: the ratio of dust mass to stellar mass as a function of specific star formation 
rate, with the same colour scheme as above. The "solar-neighourhood" and "dwarf-irregular galaxy" theoretical models of Calura et al. (2008) are overplotted 
as the open circles and the open bow-ties, respectively, with the symbols indicating the locations of the models at ages of 1, 3, 6, and 12 Gyr. The models are 
described in considerable detail in the main text, and the majority of //-ATLAS sources do not overlap with the locus of theoretical points. 



Galactic disc and centre (see Calura et al. for more detail), is indi- 
cated by the lower curve. As discussed in dCIO, the observed trend 
between specific dust mass and specific star formation rate can be 
explained as follows: dust is produced through stellar sources, its 
production rate closely linked to the star formation rate and rises 
steadily as the galaxy starts to build up stellar mass. As gas is con- 
sumed, the star formation rate declines and so less dust is formed. 
Given the destruction of dust via astration, outflows and supernova 
shocks, at this stage, galaxies can no longer replenish their dust 
content through star formation and the dust mass decreases. The 
chemical evolution models trace the evolutionary history of the 
galaxy; increasing the star formation rate, the gas mass available 
to form stars and/or the amount of dust from supernovae will drive 
the models towards the upper right in FigureQT| i.e. towards those 
//-ATLAS galaxies with the highest specific dust masses. 

The dust mass evolution as traced by the Calura/Dwek chemi- 
cal evolution models are likely to be a best case scenario since their 



assumed condensation of dust required from the heavy elements 
ejected by LIMS during their stellar wind phases and/or massive 
star supernovae are rather optimistic compared to the dust masses 
observed for stellar sources (e.g. Morgan & Edmunds 2003). In- 
deed, Dunne et al. (201 1) find it extremely difficult to explain those 
galaxies with the highest dust masses in the //-ATLAS sample 
without grain growth in the ISM as the main contributor to the 
interstellar dust budget, or a top-heavy initial mass function (see 
also Gomez et al., in prep). 

We now investigate whether the dust and star formation prop- 
erties of our sample vary as a function of stellar mass. In Figure 
[12] we show the variation of A/d us t /Altars, Md us t/SFR, and / M 
as a function of specific star formation rate binned by stellar mass, 
such that each bin has approximately equal numbers of galaxies. 
In the upper and right panels, we plot histograms for each of the 
bins in stellar mass, with the median values overlaid as vertical/ 
horizontal lines. The lower stellar mass galaxies have higher spe- 
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Figure 12. The variation of the dust to stellar mass ratio, dust mass per unit star formation rate and the / M parameter as a function of specific star formation 
rate derived from our SED fitting. The galaxies are additionally binned into three populations based on their stellar mass, with the low-, intermediate- and 
high-mass galaxies (with limits of log 10 (M/ Mq) < f0.2, f0.2 — 10.6, and > 10.6), shown by the blue, green and red squares, respectively. To the right 
of each panel we also show histograms corresponding to that mass bin's distribution relative to the vertical axis, and in the very top panel we show histograms 
of the specific star formation rate in each mass bin. 



cific dust mass, a similar trend is also observed in stacking analyses 
of optically-selected galaxies in the H- ATLAS SPIRE data (Bourne 
et al., 2012) and also in volume limited samples of very local galax- 
ies (e.g. Cortese et al. 2012), suggesting that this effect is not the 
result of selection bias in our data. 

It is also clear that the lower-mass galaxies are considerably 
more actively star-forming than the high-mass galaxies, consistent 
with the idea that the massive galaxies have consumed more of their 
available baryonic fuel through either the process of star formation 
or the aftermath of AGN feedback than their less massive neigh- 
bours (e.g. Cowie et al. 1996, Bundy et al. 2006, 2009, Hopkins et 
al., 2007, Pozzetti et al. 2010). 

The lower panel of Figure [T2] suggests that lower mass galax- 



ies have smaller contributions to their total infrared luminosity from 
dust in the ambient ISM, than their more massive counterparts 
(lower / M ). This, once more, is due to the fact that the less massive 
galaxies are undergoing proportionally more star formation and so 
the stellar birth clouds make a larger contribution to the total far 
infrared energy output. 

As in dCIO, we observe a small fraction of low-mass galaxies 
with high specific star-formation rates and high / M . The value of 
/ M in the model enters both in the optical part of the star forma- 
tion history libraries (based on the age of the stellar population and 
opacity of the birth clouds) and from the combination of the dust 
components. Clearly, there are degeneracies in the IR part of the 
SED; for example a high / M with a warm temperature for the ISM 
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component could produce a very similar looking far-infrared SED 
to a lower / M with a cool temperature for the birth cloud compo- 
nent. However, the optical colours for these two scenarios may look 
different as a result of the different stellar ages and attenuations in 
the model and so the IR data are not the only (or even the strongest) 
constraint on this parameter. 

The high / M values indicate that dust in the diffuse ISM heated 
by stars older than lOMyr dominates the far-infrared emission in 
these sources, but the high sSFR averaged over the last 0. 1 Gyr sug- 
gests that these sources have very recent star formation activity. 
If we instead use the sSFR values averaged over the last lOMyr, 
these galaxies are assigned more modest values of sSFR and shift 
to the left suggesting that the models which best fit the data are 
ones in which the star formation was recently truncated. These 
sources are all strong emission line objects, not displaying classical 
'post-starburst' spectra and so the high / M and sensitivity of model 
SSFR to the timescale of integration could indicate that they are in 
a short-lived phase transitioning from their obscured birth clouds 
to the more diffuse ISM. If global SED fitting can potentially iso- 
late sources in specific stages of evolution, this could be a powerful 
technique, however, full testing on a larger sample including a de- 
tailed analysis of the optical spectra is required in order to confirm 
this. 



4.3 Variation of SEDs within the 250 /im galaxy population 

We now analyse the shape of the averaged population SEDs as a 
function of model parameter, in order to understand the main phys- 
ical properties driving the shape of galaxy SEDs. We obtain median 
SED templates by stacking the SEDs of our galaxies according to 
their best-fit parameters, following the method described in detail 
in Appendix |Pl 

In Figure [T3] we show the median SEDs in stacks binned on 
(a) specific star formation rate, (b) dust luminosity, (c) stellar mass, 
and (d) redshift. We include only galaxies at z < 0.35 in this anal- 
ysis, i.e. we only consider the range in redshifts where we believe 
that our sample is representative. We also show the unattenuated 
starlight SED in the same bins in Figure [14] to illustrate changes in 
the underlying stellar populations fitted by the models. 

As a further check of the influence of having PACS detections 
for only a fraction of our sources, we have performed the stack- 
ing in each bin twice; once including all available data in the fit- 
ting (solid lines) and once omitting the PACS data for all sources 
(dotted lines). The number of sources in each stack is shown in 
the legend to each sub-figure, as is the fraction of these with 5a 
PACS detections. The broad similarity between the solid and dot- 
ted stacked SEDs - i.e. those including and neglecting the PACS 
data - is generally reassuring; we now discuss each set of stacked 
SEDs in detail. 



Specific star formation rate 

The most striking trends are with specific star formation rate, 
shown in Figure [T3l(a) and Figure [141(a); those galaxies with the 
highest specific star formation rates not only have the youngest 
stellar populations, but also the hottest effective dust temperatures 
manifest by the bluer optical colours, and shorter peak wavelengths 
of the FIR bump respectively. Similar results were also found by 
DCE08 for the much smaller sample of SINGS galaxies. 



Total dust luminosity 

In Figure [T3l(b) we show the relatively weak variation of the SED 
properties in our sample as a function of dust luminosity. The UV- 
optical part of the transmitted SED is similar in all but the high- 
est Ldust bin, which as shown in figure [14] is dominated by a 
younger intrinsic SED with greater reddening than the less lumi- 
nous bins. The template PAH luminosity increases markedly with 
Ldust though this area of the SED is only indirectly constrained 
by the model priors and energy balance (due to the absence of 
mid-infrared observations in this study) and so we cannot deter- 
mine how significant this is. The FIR/optical ratio increases with 
Ldust, indicating that galaxies with higher Ldust are also more ob- 
scured with a greater fraction of their bolometric luminosity being 
re-radiated by dust. The shape of the FIR peak seems to be largely 
uncorrelated with Ldust until the highest bin, when it shifts to the 
blue; thus the most luminous dust sources have warmer temper- 
atures, but it does not appear to be a monotonic trend across the 
range of Ldust probed. This change in peak wavelength is however 
consistent with the differing intrinsic (i.e. unattenuated) starlight 
SEDs shown in figure [T4l it is clear that the more dominant young 
stellar population in the most luminous bin of Ldust is associated 
with the apparently hotter dust template in the same bin. 

At this point, we need to proceed with caution in our conclu- 
sions because of the potential bias in Ldust when PACS data are 
missing (as discussed in Section 14.1.21 and Appendix[C}. The leg- 
end in figure [T3] shows the fraction of PACS-detected sources in 
each bin, and also shows in dotted lines the same stacked template 
SEDs compiled when the PACS data are neglected from the fitting 
(i.e. even for those galaxies which are detected by PACS). As figure 
1131(b) shows, the fraction of PACS-detected sources in each tem- 
plate bin of Ldust is approximately constant, and the results do not 
change when the PACS data are neglected altogether, suggest that 
they are robust to the presence/absence of PACS data. The number 
of sources in a stack in Ldust is almost always in the minority, how- 
ever, which could also be the root cause of the similarity between 
the PACS and no-PACS stacks. To check that this is not the case, 
we applied our stacking analysis to the PACS-complete sample, di- 
viding it into two luminosity bins containing approximately equal 
numbers of galaxies and stacking the best-fit SEDs derived in the 
two bins of Ldust, both including and excluding the PACS data. 
Figure [15] shows the results for the most/least luminous sources in 
the top/bottom panel, including the PACS data (in red) and neglect- 
ing them (in blue), with each stack normalised at 250 /xm. At lower 
luminosities, the lack of PACS data has no effect on the stacked FIR 
SED shape and so for bins below log 10 (Ldust / Lq) < 10.5 we can 
be confident that the lack of trend of SED morphology with Ldustis 
robust. For the bin with log 10 (Ld us t/ Lq) > 10.5, there is more 
luminosity in the mid-IR and PAH component when PACS data are 
included compared to when they are not, but this is an area of the 
SED which we cannot confidenctly discuss with the present data 
set. The FIR peak is broadened when the PACS data are included, 
but not shifted significantly, while the optical/far-infrared ratios are 
consistent, and the median optical templates are almost identical. 
We do observe an increase in the range of values (i.e. the bounds of 
the 16 th and 84 th percentiles) seen at optical wavelengths derived 
in the absence of PACS data, though they are broadly consistent. 
Indeed, these stacked SEDs are compiled from smaller numbers of 
input galaxies than any of the individual stacked SEDs in figure [T31 
and so the stacked median SEDs in figure[l5](and their percentiles, 
dotted) in particular are more susceptible to the influence of small 
numbers at all wavelengths than those in our wider study. 
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Stellar Mass 

In Figure[T3](c) we show the variation of the galaxies in our sample 
as a function of their stellar mass. Unsurprisingly, the most mas- 
sive sub-sample has a considerably more dominant old stellar pop- 
ulation in the optical wavelengths than the least massive subset, 
consistent with the unattenuated stellar SEDs stacked in the same 
way in figure [141(c). The lowest mass galaxies have broader FIR 
peaks, suggesting they have a larger warm dust component than the 
larger stellar mass sources. They also have stronger PAH emission 
in the templates, but again we caution that this is not constrained 
by data for this sample. This shift in the FIR SED shape with mass 
is likely due to the lower mass galaxies having the highest sSFR 
(i.e. "downsizing", as previously discussed) and consequently their 
average SED shows more of a warm component heated by the on- 
going star formation. 

Redshift 

Finally, we note the weak variation in the far-infrared SEDs of these 
galaxies when binned by redshift, with the panchromatic stacks 
shown in Figure [T3](d) and the unattenuated stellar componenent 
in figure[l4](d). The far-infrared SEDs of these populations all have 
similar temperatures, with the main difference between them being 
among their optical colours, with the higher redshift stacks appear- 
ing redder; figure [l4](d) shows that the intrinsic stellar populations 
in each bin are very similar, suggesting that this is due to increas- 
ing dust opacity in the UV-optical with increasing redshift. This has 
also been noted by Dunne et al. (201 1) for //-ATLAS galaxies, and 
we note that the fraction of energy emerging in the optical/FIR is 
also changing with redshift, again the higher redshift sources have 
more of their total bolometric output emerging in the FIR compared 
to those at lower redshifts, consistent with the aforementioned in- 
creased dust opacity at higher redshifts. 

The range of SEDs in each stack 

In addition to calculating the median stacked SED, we are also able 
to quantify the spread of SEDs within each bin. In Figure [16] we 
demonstrate this for galaxies binned by specific star formation rate. 
Once more, each bin has been normalised to the mean of each SED 
between 0.2 and 500 /im, and the offsets between models are ar- 
bitrary for ease of comparison. As in Figure [J_3] we present the 
transmitted galaxy templates as the solid lines, with colours corre- 
sponding to the bins of sSFR. The shaded grey regions show the 
region bounded by the 16 th and 84 th percentiles of the ensemble 
of normalised model SEDs at a given wavelength, <r(A) (see Ap- 
pendix |D]for more details); this is distinct from the much smaller 
error on the median SED, which could be used for e.g. popula- 
tion studies, or comparisons with other templates (section[5}- How- 
ever, this range in percentiles enables the reader to see the range 
of SEDs which are included in each stack, and is representative 
of how well any individual galaxy within a particular bin may be 
expected to conform to the median template. The large dispersion 
among the models at optical/ultraviolet wavelengths is due to the 
varying degrees of dust attenuation affecting the different intrinsic 
stellar populations, though as the unattenuated SEDs in figure [741 
show, the effective age of the stellar population also plays a role. 
In contrast, the comparable dispersion in the sub-millimetre wave- 
length regime (A > 80 ^tm) is due to the varying dust properties 
(e.g. temperatures, relative weights of components) of the best-fit 
models. 
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Figure 15. Stacked SEDs in two bins of Ldust derived using the PACS- 
complete sample. The top panel shows the stacks of those galaxies with high 
luminosity (log 1() (L<j U st / Lq) > 10.5, while the bottom panel shows the 
stack for the lower Ldust galaxies. The boundary values are chosen so as to 
have approximately equal numbers of galaxies in each bin, with the galax- 
ies assigned to one bin or the other according to the Ldust estimate de- 
rived using the same set of data (i.e. including/neglecting the PACS data). 
In each panel we show the stacked SEDs (solid) and 1 a scatter within the 
bin (dotted) from the best fits derived including or omitting the PACS data 
(red and blue, respectively). The far-infrared luminosities, effective temper- 
atures, and optical properties of the low luminosity stack in Ld U st(i-e. in- 
cluding/omitting PACS ) are almost identical, while for the higher Ldust hin 
there is more mid-IR and PAH luminosity when the PACS data are included, 
and the stacked SED shows a slightly broader, slightly warmer FIR peak. 

The mid-infrared wavelength range (between ~ 6 and 60 ^im), 
shows the largest dispersion, and this is expected since we lack 
good observational data at this time. Constraints in this region come 
from the upper limits at 12, 25 and 60 fj,m from IRAS (though these 
are often weak constraints), the priors fed into the stochastic far 
infrared template library (which are based on observations with 
IRAS, ISO and Spitzer), and the energy balance criterion discussed 
in Section [3TI Mid-infrared data from the WISE satellite will dra- 
matically improve this situation, especially at z < 0.2 due to the 
large fraction of 250 fim sources detected by WISE (Bond et al., 
2012). 



5 COMPARISON WITH EXISTING MODELS 

We now compare our binned SEDs to widely-used panchromatic 
SED templates, such as the empirical templates described in Chary 

6 Elbaz (2001, hereafter CE01), Dale & Helou (2002, hereafter 
DH02), or Rieke et al. (2009, hereafter R09). Each of these models 
has a strong link between the dust luminosity, and the shape of the 
SED; a link we have found to be weaker in our sample (Figure \13[. 

The CE01 templates are derived as a function of their infra- 
red luminosity, and are designed to reproduce the SEDs of existing 
WAS-selected galaxies. The data-sets used are sparsely-sampled 
and heterogeneous, using up to one hundred local galaxies at any 
given wavelength from 0.44 /im through the Infrared Space Ob- 
servatory (ISO) and IRAS bands, and out to 850/xm with SCUBA, 
albeit with no coverage between 170 and 850 /im. 
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Figure 13. The variation in the SEDs of 250 fim selected galaxies in bins of their (a) specific star formation rates, (b) dust luminosity (c) stellar mass and (d) 
redshift. The range of values for each bin, and the number of galaxies in each bin is as described in the legend for that particular plot. These median galaxy 
SEDs have been normalised to their mean flux between 0.2 and 500 fim, and offset from one another so as to reveal the different properties of each median 
SED. Whilst the vertical position of each galaxy is arbitrary, the morphlogy of each median SED is not. The differing specific star formation rates (with 
increasing specific star formation rates from the bottom of the plot upwards) are apparent in the dust properties of these median spectra in the far-infrared 
wavelengths; the most vigorously star-forming galaxies have warmer dust temperatures due to the extra source of heating in the form of a recent yield of OB 
stars. The specific star formation rate appears to be the key driver of the properties of a 250 fim selected galaxy's SED, whereas the variation between different 
bins of other properties is weaker. The stacked SEDs constructed omitting the PACS data are shown in the dotted lines, whilst those derived using all available 
data are shown as the solid lines. 



The R09 templates are based on a variety of input imaging and 
spectroscopy, including full optical photometry for 11 luminous 
and ultra-luminous infrared galaxies (LIRGs and ULIRGs) from 
the NASA extragalactic database (NED), 2MASS, IRAS, Spitzer 
and ISO, as well as the GALEXEV models from Bruzual & Char- 
lot (2003). The R09 models are binned in luminosity between 
9.75 < log 

10 (-^dust 

/L Q ) < 13.0. 

The DH02 models build upon the models of Dale et al. (2001), 
using the sample of 69 "normal" galaxies, defined according to 
their optical luminosities and Hubble types in Dale et al. (2000). 
The models have a wavelength range from 3 fim extending to ra- 
dio wavelengths, derived using data from the ISO (between 52 
and 170 fim) and SCUBA (at 450 and 850 fim) to extend the ob- 
servational constraints from 3 to 850 fim. These models represent 
"global" spectra from superpositions of local galaxy SEDs, assum- 
ing a power-law distribution for dust mass over intensity of the 
interstellar radiation field (ISRF), U, which is normalised such that 
U = 1 for the local ISRF (with values spanning 0.3 < U < 10 5 ). 
The templates are constructed such that dMd(U) oc U~ a dU, 
where Md(U) is the dust mass heated by an ISRF with intensity 



U, and the exponent (a) defines the relative contributions of the in- 
dividual local galaxy SEDs to each model spectrum. DH02 specify 
that it is those models with 1 < a < 2.5 which describe the range 
of normal galaxies. 

In Figure [T7] we overlay the CE01 models (colours) on the 
median transmitted templates from our sample (black lines with 
grey shaded regions to indicate the uncertainty on the median tem- 
plate SED following the method of Gott et al. 2001 - as opposed to 
the variation across 16 th -84 th percentiles of the stacked ensem- 
ble of galaxy SEDs as a function of wavelength in each bin as 
in figure 1 1 6b . with each overlapping set of SEDs corresponding 
to the same range in far-infrared luminosities. We normalise the 
SEDs with corresponding values of luminosity to the mean of each 
SED between 6.0 and 500 fim in wavelength, and arbitrarily offset 
them from one another in the vertical direction for ease of compar- 
ison. Whilst there is generally good agreement between the models 
at optical wavelengths, the models differ considerably in the far- 
infrared, where each CE01 model peaks at shorter wavelengths (in- 
dicative of a hotter effective dust temperature) than the correspond- 
ing /f-ATLAS stacked SED. Whilst there is some small bias in our 
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Figure 14. Stacked best-fit unattenuated (i.e. intrinsic) stellar SEDs, derived in the same z < 0.35 bins of (a) sSFR, (b) Ldust ( c ) A^stars, and (d) redshift as 
for the full panchromatic SEDs in figures[l3] The colours are identical to figures[l3](a)-(d), and each stack has been derived including all available data; once 
more, the normalisation of each SED is arbitrary, but the morphology is not. 



stacked SED in the highest two luminosity bins (see Section POl l. 
this is not an issue for the bins at log 10 (Z/dust/£©) < 10.5 and 
yet the differences between the CE01 and H- ATLAS templates 
persist. The small bias in our stacking at higher luminosities (fig- 
ure Q3) is also not large enough to account for the differences in 
the 10.5 < log 10 (i d ust/-C'o) < 11-0 bin. We also note that the 
sub-mm portion of the if- ATLAS SEDs is well determined for all 
sources due to the high quality Herschel SPIRE data that form the 
foundation of this study; the differences at these wavelengths rela- 
tive to the CE01 templates persist across the full luminosity range 
of our sample. 

The largest disagreement between the two sets of models is 
at mid-infrared wavelengths, however as lack observations in this 
region of the SED the current mid-IR discrepancy with other tem- 
plates is not significant; a full analysis with WISE data will be re- 
quired to see if these differences persist. 

In Figure [T8j we compare our templates with those of R09, 
this time normalised between 100 & 500 /im. The R09 templates 
are supplied only at A > 4 /im, and so the range in wavelength val- 
ues is smaller than for the other three sets of models. We find that 
the R09 models have similar dust temperatures to our models, as 
demonstrated by the similar peak wavelengths of the far-IR SED, 
but that they all have considerably brighter mid-infrared emission 
than our templates. Again, we cannot comment further on this dis- 
crepancy until we have been able to consider the WISE data in the 
fitting. 



In Figure [79l we compare our models with the DH02 tem- 
plates, which are empirically constrained at A > 3 /im. Rather 
than limit the comparison to the range of values of a which DH02 
suggest span the range of normal galaxies in their input sample 
(1.0 < a < 2.5, magenta in Figure [79t, we also compare our 
SEDs with the more quiescent range of DH02 templates spanning 
2.5 < a < 4.0 (light blue in Figure fT9l. We normalise the models 
at 250 /im for ease of comparison. Whilst there are DH02 models 
which can match our templates at the highest dust luminosities, at 
lower values the DH02 models suggest the presence of consider- 
able hot dust components that are not observed in our sample of 
250 /jm-selected galaxies, and the a > 2.5 models are required to 
reproduce the far-IR temperatures that we observe. Since the DH02 
models are binned by a rather than dust luminosity, we indicate the 
luminosity of our models being compared using the colours as in- 
dicated in the caption to Figure[T9] with the DH02 models overlaid. 

These templates may be useful for studies of FIR selected 
samples in the z < 0.5 universe. They are also useful ingredients 
for any evolutionary model which purports to explain the low red- 
shift FIR/sub-mm populations detected by Herschel. They do not, 
however, appear to be representative of the high redshift (z > 1) 
population detected in W-ATLAS (e.g. Lapi et al., 201 1). The tem- 
plate chosen for any SED fitting should always match the sample 
under investigation as selection wavelength and redshift can have 
an important impact on the SED types prevalent in a sample. 

Irrespective of the templates to which we compare our model 
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Figure 16. In addition to Figure [T3l it is possible to show the variety of different SED models within a bin of e.g. specific star formation rate. The median best 
fit transmitted SEDs are shown by the coloured solid lines corresponding to the bins of sSFR, with the la spread of SEDs going into the stack shown by the 
shaded areas bounded by the black lines. In making these SEDs, each component SED in a particular bin of specific star formation rate has been normalised 
to its mean between 0.2 and 500 /im and the resulting median SEDs (plotted) have been offset from one another to prevent them from overlapping. The large 
variety of best fit models at mid infrared wavelengths (between ~ 6 — 60 fim) is due to the lack of observational constraints in this part of the spectrum. 



SEDs, the comparative lack of hot dust in our stacks is striking. 
Though we plan to investigate these details in future using WISE 
observations, the contribution of the mid-infrared to the total dust 
energy budget is not dominant; the impact of the WISE data near 
the peak of the far-infrared SED is likely to be small. 

We intend to make these new template SEDs, binned accord- 
ing to their properties available to the community for further anal- 
ysis, and application to other data sets via the //-ATLAS websittQ 
and the author's webspac^E 



6 CONCLUSIONS 

We have determined SEDs for a total of 1402 250/im-selected 
galaxies from the Herschel-ATLAS science demonstation cata- 
logue with reliable counterparts and matched aperture photometry 
from the it- to K -bands from the GAMA database. We also include 
far and near ultraviolet data from the GALEX-GAMA survey, as 
well as the //-ATLAS data from PACS and SPIRE. Of these 1402 
galaxies, 1289 are well described by the model of DCE08, and we 
use these SEDs and the model parameter probability density func- 
tions derived from the energy balance SED-fitting, to determine the 
properties of these 250 fim selected galaxies out to z = 0.5. 



http : / / www . h— atlas ■ org| 



http : / /star . hert s . ac . uk/ ~dsmith/| 



(i) Studies of the colours of galaxies in our sample, and a suite of 
simulations, suggest that our sample is representative of the broader 
population of 250 galaxies out to z < 0.35. 

(ii) The average //-ATLAS galaxy in our sample has a star for- 
mation rate of ~4.0 M yr _1 , L du3t « 6.4 x 10 10 L , and a dust 
to stellar mass ratio of ~0.4 per cent, while the median redshift is 
z = 0.24. 

(iii) Our results support the idea that IRAS and //-ATLAS se- 
lected galaxies in the local Universe are different populations. Due 
to its lack of sensitiity and short selection wavelength, IRAS prefer- 
entially selected galaxies with larger warm dust content, and con- 
sequently these galaxies are more luminous in the infrared for a 
given mass of dust. The //-ATLAS selection at 250 /im is less bi- 
ased towards strongly star forming objects over the same redshift 
range because of the longer selection wavelength and far superior 
sensitivity compared to IRAS. IRAS misses a population of massive 
dusty galaxies with colder dust temperatures, as was shown previ- 
ously by Vlahakis, Dunne & Eales (2005). 

(iv) The correlation between star formation rate and dust mass 
presented in da Cunha et al. (2010) is also present in this sample, 
although Herschel ATLAS-selected galaxies contain larger dust 
masses for a given star formation rate compared to the IRAS se- 
lected sample of dClO. There is also a correlation between specific 
dust mass (A/d us t/A/st ars ) and SSFR, which is not well reproduced 
by simple chemical and dust evolution models. 

(v) The specific star formation rate of lower mass galaxies 
(log 10 M B tar/Mo < 10.2) is higher than that of the most mas- 
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Figure 17. Comparison of our median, transmitted, and luminosity-binned 
templates (solid black lines, with the uncertainty on the median template 
SED in each stack - defined using the median statistics method of Gott et 
al., 2001 - indicated by the shaded regions) with the templates of Chary & 
Elbaz (2001 - CE01, shown in colour and with the luminosity bounds in- 
dicated by the legend). Each set of transmitted models with the same range 
in dust luminosity (indicated by the colours) has been normalised to have 
the same mean between 6 & 500 fim, and artificially offset in the vertical 
direction to prevent the models from overlapping. 
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Figure 18. Comparison between the Rieke et al. (2009) and our template 
SEDs binned according to luminosity. As in Figure [171 each set of SEDs is 
normalised between 100 & 500 fim (the Rieke et al models are not defined 
in the optical), and arbitrarily offset, with our templates and uncertainties 
indicated by the solid line and dashed areas, respectively. Once more the lu- 
minosities of the R09 templates are shown by the colour, tabulated accord- 
ing to the legend. Although the R09 templates are reasonably similar to the 
new //-ATLAS templates at the lower dust luminosities, at higher dust lu- 
minosities they become increasingly warmer than the //-ATLAS templates. 



sive galaxies in our sample (those with log 10 A/ Star / Mq > 10.6) 
at all redshifts, supporting previous results that lower mass galaxies 
dominate the star formation rate density in the local universe. 

(vi) Stacks of SEDs show that sSFR is the strongest galaxy prop- 
erty driving the SED shape across both the UV/optical and FIR, as 
first noticed in the smaller sample of DCE08. Trends with idust 
are much weaker since smaller mass galaxies will have low Ld us t 
and yet could have the highest values of sSFR. We see a signficant 
trend in this sample for galaxies to have more obscured optical/UV 
SEDs and higher reprocessed fractions with increasing redshift. 

(vii) Existing templates for panchromatic SEDs of galaxies 
show shorter FIR peaks and excess mid-IR emission compared to 
median stacked SEDs of galaxies in our //-ATLAS sample (binned 
by idust) although the mid-IR discrepancy is not significant at this 
time due to our lack of mid-IR data to constrain this part of the 
SED. Templates from Rieke et al. (2009) are the closest match to 
ours in terms of the FIR properties although they still predict a 
warmer FIR peak at the highest luminosities compared to our find- 
ings. We provide a new set of panchromatic SED templates from 
the UV-sub-mm, to enable more representative studies of dusty 
galaxies in the local Universe in the Herschel era. 

(viii) Data from the WISE satellite, which covers the wavelength 
range between 3 and 23 /J,m will provide valuable constraints to the 
mid-IR and PAH features, as well as the hot dust component of 
these local galaxies. It will be interesting to see if the differences 
between templates in the mid-IR region persists when these data 
are included in the fitting. 
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Figure 19. A comparison between our templates and those of Dale & Helou 
(2002). The DH02 templates are shown in magenta for the range of a val- 
ues thought by DH02 to describe normal galaxies, and in light blue for the 
most quiescent galaxies in the DH02 library (with 2.5 < a < 4.0). The 
luminosities of our templates are shown by the colour of the solid lines 
bounded by the shaded regions indicating the uncertainty associated with 
each template. For ease of comparison, each set of models being compared 
is normalised around 250 (im, and arbitrarily offset in the vertical direction. 
It is clear that the DH02 templates contain considerably hotter dust than we 
observe in 250 fim selected galaxies in //-ATLAS. 
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APPENDIX A: STACKING PDFS AND CHOICE OF COLD 
TEMPERATURE PRIOR 

Al Derivation and interpretation of stacked PDFs 

In this paper we make considerable use of stacked PDFs, which 
are our best estimates of the distribution of the values of a given 
parameter amongst the sources in a sample, convolved with our 
ability to constrain them. To see why this latter point might be im- 
portant, imagine we have individual PDFs for a certain parameter 
which has very weak constraints, all the PDFs will therefore ap- 
pear flat. The distribution of the medians of these PDFs (i.e. the 
median-likelihood estimates for that particular parameter) will be 
centered near the middle of the range (since the PDFs are all in- 
dividually flat) and will have a narrow scatter (since all medians 
are almost the same). If we only considered the median-likelihood 
values, we might naively assume that we know the parameter distri- 
bution for the population quite accurately, and that there was little 
scatter within the population, even though in truth we merely had 
little ability to constrain that parameter. This is an extreme example 
and does not apply to the parameters we are exploring in this paper, 
but it illustrates why we wish to present the stacked marginalised 
PDFs and not simply the distributions of the median likelihood val- 
ues. 

To derive values of stacked PDFs, we start with the values 
in each bin of the ensemble of PDFs that we wish to stack. For 
the value of the stacked PDF in each bin we use the mean of the 
ensemble of values in that bin. To estimate the error associated with 
the derived stacked PDF in each bin, we simply use a symmetric 
value corresponding to the mean of the 16-84th percentiles of the 
cumulative frequency distribution of the values for each galaxy in 
that bin. 

A2 Choice of T co i d prior distribution 

When choosing a prior distribution for a particular parameter in 
the stochastic libraries of SEDs, we must ensure that our choice 
of prior does not bias our results. This was of particular concern 
for the prior on T co id since the stacked PDFs do not always show 
a peak in the range of the prior, but for some samples, increase 
towards the lowest bounds of the prior. We wanted to determine 
whether or not this was an indication that we should use a broader 
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Figure Al. The stacked PDFs (solid line, with error bars) for those 59 
galaxies with PACS 160 /^m detections with 5250 > 120 mJy at z < 0.20. 
We overplot the median values of the individual PDFs for these galaxies 
(dotted histogram), along with the best-fit Gaussian approximation of the 
PDF. The deconvolved version (described in the text) is shown in the dot- 
dashed line, and has a standard deviation of 3.54 ± 0.41 K. 



prior on T co id. To address this issue, we considered the properties 
of the galaxies in our PACS-complete sample (see Section 14.1.21 . 
which as we showed in section l4~, 1.2l and figure[6] are representative 
of the full range of colours in //-ATLAS sources, and by virtue of 
their being detected in our PACS data, have our best constraints on 

Tcoid- 

We stacked the PDFs for these galaxies, and the results are 
shown in Figure [ATI in which the best fit Gaussian approximation 
to the stacked PDF is shown as the dashed line, and a histogram of 
the median values of each individual PDF that went in to the stack 
in dotted lines. The best-fit Gaussian model of the stacked PDF 
has a standard deviation of 3.84 ± 0.41 K. In order to assess the 
true range of T co id that is present in our sample, we determined the 
deconvolved best-fit Gaussian, by subtracting the mean of the indi- 
vidual la errors on T co id for each galaxy (1.47 K±0.19) in quadra- 
ture, leaving a la uncertainty on the range of T co id in our sample as 
3.54 ± 0.41 K. The best fit deconvolved Gaussian is shown as the 
dot-dashed line in Figure [ATI We find that 85 percent of the true 
T co id PDF lies within the bounds of our temperature prior, with 
only ~six per cent of the PDF colder than 15 K, and approximately 
9 per cent warmer than 25 K. 

In that case, why not make a T co id prior which is wider and 
therefore encompasses the full range of temperatures possible in 
the cold ISM? The reason not to do so is because of the strong non- 
linear dependence of the Md us t parameter on T co id when the value 
of Tcoid is below ~ 15 K. At such cold temperatures, the SPIRE 
bands no longer sample the Rayleigh- Jeans part of the SED (where 
A/dust scales relatively linearly with T co ki ) but are nearer the peak 
where the dependence on temperature is steeper. 

The fitting becomes more prone to errors on the determination 
of Tcoid at low values, since the energy balance is not much affected 
by the addition of very cold dust to the SED (which contributes lit- 
tle to the total Ldust). Dust colder than 15 K is essentially invisible 
to our model (or any other for that matter), which combined with 
the added sensitivity of mass to temperature at low T co id results in 
an asymmetry in the error on the dust mass; we have a larger over- 
estimation of the mass when T co id is underestimated (statistically 



© 2012 RAS, MNRAS 000.iTll26l 



Herschel ATLAS: SEDs 23 



likely to happen 50 per cent of the time) compared to the size of 
our underestimate in Md us t when T co id is overestimated. 

To demonstrate this, we created a hypothetical library of 1000 
model galaxies with Gaussian temperature and mass distributions, 
and simply re-estimated the dust mass after adding on simulated 
measurement errors to the true temperature distribution, using the 
relationship between temperature and dust mass for a modified 
greybody emission profile shown in Equation 3 of Dunne et al. 
(201 1). In these simple simulations, measurement errors introduced 
a systematic anti-correlation between the estimated temperature 
and estimated dust mass. In Figure lA2l we show a histogram of 
the difference in our median estimated dust masses, derived us- 
ing the broader and narrower T co id prior distributions. The effect is 
stronger at the coldest temperatures, which occur more often with 
the broader T co id prior, with some dust masses being overestimated 
by in excess of 0.5 dex. In order to limit the effects of this unphysi- 
cal bias towards large dust mass estimates in the colder galaxies in 
our sample, we decided to use the narrower prior on Tcoid- 

We confirm the effects on T co id by re-running the fitting using 
a wider T co id prior (10-30 K) on sub-samples of galaxies which 
represent typical selections within the main analysis and compare 
these results to those using the T co id prior from DCE08 (15-25K). 

We split our sample into five subsets, limited to z < 0.2 to 
limit the possible influence of cosmic evolution on our results: 

• galaxies in our PACS-complete sub-sample 

• galaxies detected at > 5a in both PACS bands and > 3a at 
350 & 500 /im; 

• galaxies with at least one PACS ^ 5a detection; 

• galaxies detected at ^ 3a at 350 /im; 

• all galaxies in our sample. 

Of course, all galaxies discussed here are detected at ^ 5a at 
250 /im, since this is how our sample is defined. 

In Table IaTI we show the median-likelihood values of T co id 
derived from the PDFs for each sub-sample. None of the stacked 
PDF median-likelihood T co id values varies by more than ~0.6K 
when the broader prior is used rather than the narrower prior; the 
impact on our global estimates of T co id is therefore minimal. When 
we use the broad prior on T co id, we find that 13 per cent of our 
sample have T co id values of < 15 K, approximately in line with 
the 6 per cent that we expect from studying the PACS-complete 
sample. 

When we consider the median temperature estimates with 
each prior in table lAll we see that the second two samples - which 
require 5a PACS detections - have higher median values due to the 
removal of the coldest sources from the sub-sample; these sources 
are generally undetected in our comparatively shallow PACS data. 
The last two sub-samples in table lAll have lower values of T co id, 
and are mostly undetected by PACS; in particular, the 350 se- 
lection criterion preferentially picks the colder sources in the cata- 
logue, with SEDs peaking at longer wavelengths. 

In summary, we have used the DCE08 15 K < T co id < 25 K 
temperature priors for the following reasons: 

• Simulations and studies of the PACS-complete sample suggest 
that we can explain our data adequately using the DCE08 prior, 

• The narrower prior also limits the impact of T co id errors on 

Mdust, 

• Detailed studies using a variety of multiwavelength data and 
modelling techniques see no evidence for substantial cold dust 
components with temperatures lower than 15 K in the kind of dust 
rich spiral galaxies being investigated here (e.g. Dunne & Eales 
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Figure A2. The difference in the dust mass estimates that we determine 
for the galaxies in our sample depending on whether we use the broader 
(10 < T co id < 30) or narrower (15 < T co id < 25) prior. The asymmetry 
in this histogram is caused by a greater sensitivity to temperature at cold 
temperatures. 



2001, Popescu et al. 2002, Vlahakis, Dunne & Eales 2005, Draine 
et al. 2007, Willmer et al., 2009, Bendo et al, 2010, Boselli et al. 
2010, Kramer et al. 2010, Bernard et al. 2010). 



APPENDIX B: DEFINING THE GOOD FITS 

We took a selection of best-fit model SEDs, and varied their pho- 
tometry according to a set of Gaussian distributions with a median 
of zero and a standard deviation equal to the minimum photometric 
error in each band (i.e. the values added in quadrature, as defined 
in Section|2j- These values were chosen since they were the domi- 
nant source of error for bright sources. We then removed a sub-set 
of photometry to reflect the heterogeneity within our real data set, 
and re-calculated their properties 1000 times each. The resulting 
histograms of \ 2 values (e.g. Figure IBTt enabled us to estimate the 
number of degrees of freedom given that particular sub-set of pho- 
tometry, by performing a simple minimisation of Equation IB II to 
the derived \ 2 probability density functions: 



P oc 



1 



2("dof/2)r(JV do f/2) 



(x 2 ) 



2\(iV dof /2)- 



1 -X 2 /2 

exp A ' 



(Bl) 



where Y{x) represents the Gamma function, JVdof the number of 
degrees of freedom, and \ 2 the median \ 2 value for each bin. Due 
to the relatively small redshift range covered by our sample, we 
average each solution, and determine that the relationship between 
the number of degrees of freedom (JVdof) and the number of pho- 
tometric bands with measurements (A^ands) is given by Equation 



(B2) 



A dof « (-2.820 ± 0.745) + (0.661 ± 0.132)/V bands 

+(7.91 ±5.50 x 10~ 3 )iV b 2 ands . 

With this information, we are then able to use Equation IB II in 
conjunction with the number of degrees of freedom estimate from 
Equation lB2l to determine a 99 per cent confidence interval on % 2 - 
Those galaxies outside the interval on x 2 have less than 1 per cent 
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Table Al. The effects of our choice of prior distribution on our median-likelihood estimates of T co id . The left column indicates the selection of each different 
sub-sample, in addition to the requirement that all sources have 5<r detections in the SPIRE 250 ftm band. The middle two columns indicate the T co id 
estimates determined for each sub-sample with the two prior distributions, while the right hand column shows the difference. Each sub-sample has been 
limited to z < 0.2, in order to limit the impact of evolution with redshift on these values. 
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Figure Bl. The distribution of \ 2 values for the one thousand galaxy sim- 
ulations referred to in the text. The shape of the resulting histogram is well- 
described by Equation lBll and depends on only one parameter, the number 
of degrees of freedom for that particular combination of input photometry 
(Ndol )■ In this way, N^ a f can then be determined for each combination of 
photometry by simply minimising Equation IB II By fitting a quadratic to 
these results, we may relate A^of and A Dan ds- The resulting relationship 
is shown in Equation lB2l The combination of Equations IB 1 l & lB2l allows us 
to define 99% confidence intervals and in turn determine which galaxies in 
our sample are well-described by our model. 



chance of being consistent with our model. In this way we may re- 
move those "bad fits" from further analysis. We derived equation 
IB2l using galaxies with between 6 and 19 detections, so these val- 
ues constitute the bounds on iVbands over which we believe it is 
valid. 

Finally, in comparing our simulations to our real sample, we 
note that in contrast to our simulations, the photometric errors used 
in our SED fitting are not strictly Gaussian due to problems as- 
sociated with e.g. deblending or calibration issues, particularly in 
regions neighbouring saturated stars, however the difference is not 
expected to be large. 



APPENDIX C: 
FITTING 



ADDITIONAL TESTS FOR BIAS IN THE 



Due to the complexity of the energy-balance SED fitting method, 
we wanted to perform additional checks for bias in the derived pa- 
rameters, and determine which parameters depend on one another. 
In turn, we calculate the variation on our key parameters 



(Ldust, Mdust, fit, Altars, SFR and sSFR) introduced in our 
"PACS-complete" sample when the PACS data are included and 
when they are omitted from the fitting, as a function of each of 
the derived parameters included in this analysis. We consider the 
difference in each parameter, Q (e.g. Ldust, sSFR, etc): 



AQ = Q 



(CI) 



where Q PACS and Q noPACS are the median likelihood estimates 
for a particular source including or ignoring the PACS information 
in the fitting. In this way we are able to probe for a skew in the bias 
between estimates of a parameter e.g. do we over-estimate Ldust 
for particular Ld us t when galaxies aren't detected by PACS? Even 
though we have used our stacked PDFs to probe for overall bias 
and found that it is minimal for all of the key parameters that we 
discuss in this paper (they are detailed in tabled, it is still possible 
that AQ is skewed for some parameter combinations. In figure ICTl 
we plot these values along with their error bars. The individual data 
points are shown as red circles, with error bars in the horizontal 
direction derived in the absence of the PACS data. The blue solid 
line indicates zero offset, while the green solid line indicates the 
median of the stacked PDF for all galaxies in the "PACS-complete" 
sample that we explore here. 

In figure |CT| (a), we consider possible bias in Ldust- We find 
that there is a skewed bias between ALdust and ALju.rt ACS ; though 
the scatter is large, these tests suggest that Ld us t is underestimated 
at high Mdust and underestimated at lower values in the absence 
of PACS data, though the offset averaged across the whole range 
of Mfaft is small. ALdust appears broadly unbiased across all 
values of Ld u PACS , while there is little evidence for any bias with 
respect to M sta rs, fit, (specific) SFR or redshift. There is evidence 
for bias in L du st with respect to Mlust/Ldust and ALjust/M star , 
though these biases are dominated by the bias in Afdust already 
discussed. 

Figure ICTl (b) shows the variation in AALjust as a function of 
the same key parameters, suggesting that there may be a tendency 
to overestimate Mdust in the absence of PACS data at the highest 
dust masses probed by this study and there is a small tendency for 
dust masses to be overestimated in the absence of PACS data on 
average (comparing the blue and green horizontal lines in figure 
ICTl The related parameters M d 1 °S ACS /L5° PACS and M^°f t ACS / 
M£t£r S ) also show evidence for skew. 

Figures IC 11 (c) and (d) show the variation in A / M and 
AALtars; it is clear that they are unbiased with respect to the other 
model outputs that we use in this paper, as probed by the "PACS- 
complete" sample. Though the spread in ASFR and AsSFR is 
larger, reflected in the larger error bars in figure ICTI fe) & (f), it is 
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(e) SFR (f) sSFR 

Figure CI. Variation in (a) Ldust. 0>) Afd u8 t, (c) / M , (d) M a tars, (e) SFR and (f) sSFR as a function of each of the parameters discussed in this paper for 
the PACS-complete sample described in Section l4. 1.21 The data points are shown as red circles, with horizontal error bars indicating the uncertainties derived 
neglecting the PACS data, and the vertical position derived according to equation lCll considering the asymmetric errors on the median likelihood values for 
each parameter. The blue horizontal line indicates AQ = whilst the green horizontal line denotes the difference between the median values of the stacked 
PDFs for the whole PACS complete sample when the PACS data are included and when they are omitted. 
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Figure Dl. Individual best-fit SEDs (grey points) are first normalised be- 
tween 0.2 and 500 fj,m, and then the median (the red line) and 16 th and 84 th 
percentiles of the ensemble (blue lines) are calculated as a function of wave- 
length. This particular example shows the stacked SEDs of those 216 galax- 
ies with good best-fit models at z < 0.35 and with 10.0 < log 10 (iduat/ 
-^solar) < 10.5. 



difficult to discern any skew or bias between these parameters and 
the other model outputs. 



APPENDIX D: STACKING SAMPLES OF SEDS 

In order to calculate the median of an ensemble of SEDs, we 
first normalise each individual SED to the mean between 0.2 and 
500 (im (in units of \F\), such that each template is given equal 
weighting in the stack, then take the median of the ensemble of 
normalised SED values in each wavelength bin. We also determine 
the 16 th and 84 th percentiles of the cumulative SED-distribution 
as a function of wavelength. These values provide a measure of the 
spread in the SEDs of the galaxies which go into the stack, in con- 
trast to the estimated error on the median template shown in figure 
|16| which is determined using the median statistics method of Gott 
et al. (2001). Our method of determining the range of values in our 
SED stacks (i.e. the 16-84 th percentiles) is illustrated in Figure lDTl 
in which the individual best-fit SEDs that go into the median tem- 
plate (in this case for 216 galaxies with good fits, at z < 0.35, 
and with 10.0 < log 10 (Ldust/isoiar) < 10.5) are shown in grey, 
while the median and the afore-mentioned percentiles are overlaid 
in red and blue lines, respectively. 

Note that for the purposes of calculating these stacked tem- 
plates, we bin according to the best-fit values returned from the 
fitting code rather than the medians of the PDF that we use for anal- 
ysis elsewhere. This distinction is noteworthy since we determine 
only the best-fit model SEDs for each galaxy (as opposed to the full 
PDF at each wavelength). Whilst there is generally excellent agree- 
ment between the best-fit and the median-likelihood estimates of 
any given parameter, the two may differ in individual cases, adding 
unrealistic outliers in the stacked SEDs if median-likelihood values 
are used for these purposes. 

Whilst this method of stacking SEDs determines templates de- 
signed to be representative of a typical galaxy in a given sample (or 
sub-sample) of galaxies, it is important to note that we do not ex- 
pect such templates to reproduce the total cosmic spectral energy 
distribution (e.g. Hill et al., 2010, Somerville et al., 2012, Driver et 
al., in prep), which would require calculating the sum of the emer- 
gent SEDs without prior normalisation. We leave further discussion 
of this topic for a future publication. 



© 2012 RAS, MNRAS OOO.QJTJg] 



